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The theory of orthogonal polynomials on the unit circle is developed 
for a general class of weights leading to systems of recurrence relations 
and derivatives of the polynomials and their associated functions, and to 
functional-difference equations of certain coefficient functions appearing 
in the theory. A natural formulation of the Riemann-Hilbert problem is 
presented which has as its solution the above system of orthogonal poly- 
nomials and associated functions. In particular for the case of regular 
semi-classical weights on the unit circle w{z) = rijLi(-^ ~ '^ji^)Y^ y con- 
sisting of m G Z>o singularities, difference equations with respect to the 
orthogonal polynomial degree n (Laguerre- Freud equations) and differen- 
tial equations with respect to the deformation variables Zj (t) (Schlesinger 
equations) are derived completely characterising the system. It is shown 
in the simplest non-trivial case of m = 3 that quite generally and simply 
the difference equations are equivalent to the discrete Painleve equation 
associated with the degeneration of the rational surface D^^ ^ — > Dg^' and 
no other. In a three way comparison with other methods employed on 
this problem - the Toeplitz lattice and Virasoro constraints, the isomon- 
odromic deformation of 2 x 2 linear Fuchsian differential equations, and 
the algebraic approach based upon the affine Weyl group symmetry - 
it is shown all are entirely equivalent, when reduced in order by ex- 
act summation, to the above discrete Painleve equation through explicit 
transformation formulae. The fundamental matrix integrals over the 
unitary group U{N) arising in the theory are given by the generalised 
hypergeometric function 2^1^'- From the general results flow a number 
of applications to physical models and we give the simplest, lowest order 
recurrence relations for the gap probabilities and moments of character- 
istic polynomials of the circular unitary ensemble (CUEjv) of random 
matrices and the diagonal spin-spin correlation function of the square 
lattice Ising model. 
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1. Introduction 

With —TV < < TV, zi := e*^' the unitary group U{N) with Haar (uniform) measure 
has eigenvalue probability density function (see e.g. |18l Chapter 2]) 

(1-1) (2^^ n i-'^-^A'- 

^ ' l<j<k<N 

Our interest is in averages over U{N) of class functions it;(f/) which have the factorization 
property YliLi ^i^i) for {^i, • ■ • > ^n} G Spec(L''(A'')). Introducing the Fourier components 
{wi}i=o,±i,... of the weight w{z) by w{z) = J^H-oo due to the well known identity 

N 

(1.2) (n^(-^')) = det['u;i_j]ij=i,...,]v, 

we are equivalently studying Toeplitz determinants. As an explicit example consider the 

unitary average 

(1.3) 

with X(2-ci> tt) ~ f Oi £ {tt — (t>-,i^) and X(2-c/) tt) ~ ^ otherwise. For special choices of 
the parameters 1)1. 3|l occurs in a variety of problems from mathematical physics. Thus the 
case (u;i,ciJ2,m) = (0,0,0) is the generating function for the probability that the interval 
(tt — (j), n) contains exactly k eigenvalues in Dyson's circular unitary ensemble (which is 
equivalent to the unitary group with Haar measure), while the case (coi , u!2 , fJ.) — (1,0, 1) 
is (apart from a simple factor) the generating function for the probability density function 
of the event that two eigenvalues in the circular unitary ensemble of (A'^ + 2) x {N + 2) 
matrices are an angle <j!> apart with exactly k eigenvalues in between. The case ^ = 2, 
a;2 = 0, /X = oil = 1/2 of 11.31 1 corresponds to the density matrix for the impenetrable 
Bose gas |^^]. Furthermore in the case 5 = one sees that 11.31 includes as special 
cases 

AT 

(1-4) (l[zj^'\l + z,r/'{l + k-'z,f/') , l/k'<l, 

\ fj^ I C7(JV) 



(1.5) (n(i + iM)"'(i + g'20") , a'<i- 

1=1 ^ ' 

The average 11.411 is equivalent to the Toeplitz determinant given by Onsager for the diag- 
onal spin-spin correlation in the two-dimensional Ising model |42|. while 11.511 occurs as a 
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cumulative probability density in the study of processes relating to increasing subsequences 
[7IIT1). 

In |22l the average (11.311 was characterized as a r-function in Okamoto's Hamiltonian 
theory of the sixth Painleve (Pvi ) equation, up to a change of variables and/or multipli- 
cation by an elementary function. The r-function r(t) is defined to be 

(1.6) ff^llogr, 

where the Hamiltonian H{q,p; a, i) is a rational function of the co-ordinates and momenta 
q,p and of parameters a — (ao, ai, 02, as, 04), {ao -I- Qi -I- 2a2 -1-03-1-04 — 1), and 
independent variable (deformation variable) t. The dynamics Tt : g(0),p(0) q(t),p(t) is 
governed by the Hamilton equations 

, > dq _ dH dp _ dH 

^ ' ' ~ 'dt ^ ^~d^ 

where eliminating p{t) gives the sixth Painleve equation in q{t). The above identification 
implies that the logarithmic derivative of Tn with respect to t is an auxiliary Hamiltonian 
for the Pvi system and it was shown in f22] that it satisfies a difference equation related 
to a particular discrete Painleve equation with respect to increments in unit amounts of 
one of the parameters a (or in particular fi). One of the objectives of this paper is to 
derive discrete Painleve type recurrences for Tn directly with respect to increments in A*' 
only. In the algebraic approach PO) the strategy is to choose a particular shift operator 
or Schlesinger transformation L constructed from compositions of fundamental reflection 
operators and Dynkin diagram automorphisms of the Pvi symmetry group Wa{D'^^) — 
(so, Si, S2, S3, S4, ri, r3, r^) and use it to generate sequences of parameters and dynamical 
variables L" : {a,q,p,H,T} ^ {an,qn,Pn,H„,Tn} for n = 0,1,.... This operator is 
chosen to increment the parameters such that only TV A'^ -I- 1 and thus rjv is essentially 
the average l|1.3^ . In all our applications the initial point a{N — 0) is located on a 
reflection hyperplane in the space a £ C* and this characterises the entire sequence as 
classical solutions to the Pvi system. The initial member tq — 1 and the first nontrivial 
member ti is a solution of the Gauss hypergeometric differential equation. 

The discrete Painleve equation that is fundamental in the Pvi system is that associated 
with the degeneration of the rational surface D^'' — > Dg^' 

,(/n + 1 - 02) (/n + 1 - ao - 02) 



/n(/n + as) 

Qi a4t 



(1-8) g„+ign=t- 

(1.9) /n + /n-l = -Q3 + ■ ^ , 

where ai 1-^01 + 1,02 1-^02 — 1, 04 t-^ a4 + l as n ^ n + 1. We label the particular discrete 
equations that arise according to the unambiguous algebraic-geometric classification of 
Sakai |44| by association with a degeneration of a particular rational surface into another, 
rather than the historical names employed (discrete fifth Painleve equation, dPv). In 
its full generality the average 11. 3t was first characterised in terms of 11.811 and 11.91 in 
|20| . In fact this discrete Painleve equation (but different in detail than the ones we will 
present below including that in (2()) 'l for the average H1.5|l is already known from the work 
of Borodin [12 j. In contrast recurrences for the average l|l.i-i|l with ^ = have also been 
obtained recently by Adler and van Moerbeke from their theory of the Toeplitz lattice 
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and its Virasoro algebra that do not appear to relate to and ljl.9|l Another 

of our objectives is to show that indeed the recurrences of are of the discrete Painleve 
type, by deriving explicit transformation formulae between them and ll.8ll . ll.9t . 

Because we are dealing with Toeplitz determinants with symbols w{z) satisfying cer- 
tain analytic conditions it is immediate from the theory of Szego, Geronimus and oth- 
ers that orthogonal polynomial systems on the unit circle {<;/>„ (z)}55^^o with respect to 
such weights are relevant. For example if £ [0,2tt), uji,uj2, £ R, ^ < 1, and 
2tJi > -l,2^i > -1 then w(e'*) where 



w{z) = ^-'^-"(1 + (1 + tzf" 



1 e ^ {tT - (t>,Tv) 

e G (tt -(?!., tt) 



is a real, positive weight defining a measure with an infinite number of points of increase 
and thus an orthogonal polynomial system on T exists with respect to this weight by 
Favard's theorem. To characterise the averages Ijl.^^ using orthogonal polynomial theory 
we have found it necessary to substantially develop the general theory of such systems. To 
a large extent this task has been completed for orthogonal polynomial systems defined on 
the line in the works of Bauldry 0, Bonan and Clark [ill] , Belmehdi and Ronveaux 
Magnus |88l 18711!^ but had remained incomplete for those systems on the unit circle |27| . 
To this end we have derived closed systems of differential relations for the polynomials, 
their reciprocal polynomials {0* (z)}55^oi s^nd associated functions {e„(2)}55^oi {^n{^)}n'=o 
in Proposition 12.11 In the notation of 12.521 and Corollary 2.3 let 



and set 



0n(2) e„{z)/w{z) 

0n(^:) -^*n{z)/'w{z) 



dz 



Entries in the matrix A„ are fixed by four coefficient functions Q„{z), Sl^(z), 0„{z), On{z) 
in 12.5411 and complete sets of difference and functional relations for these coefficient 
functions are given in Proposition 12.21 and Corollarv l2.2l We also formulate a 2 x 2 matrix 
Riemann-Hilbert problem in Proposition 12. 31 for general classes of weights which parallels 
the case for orthogonal polynomials on the line [2811161 [T71ll5 l and whose solution is simply 
related to Yn. For our particular applications the weight 111.311 is a member of the regular 
semi-classical class 

m 

w{z) = l[{z-z,{t)r^, ec, 

with an arbitrary number m of isolated singularities located at Zj{t). A key feature of 
such weights is that 

Id 2V{z) 
-w{z) — 



w{z) dz ^ ' W{z) ' 
where the polynomials degV^(2) < m,degW(z) — m. The coefficient functions for regular 
semi-classical weights are polynomials of z with bounded degree degf2„(2:) — degf2* (2:) — 
m — l,deg0„(2:) — deg0*(2:) — m — 2 (see Proposition 13.111 . In addition evaluations 
of these functions at the singular points satisfy bilinear relations (see Proposition 13.21 
which lead directly to one of the pair of coupled discrete Painleve equations. Deformation 
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derivatives of the linear system of differential equations above with respect to arbitrary 
trajectories of the singularities are given in Proposition 18.31 which can summarised as 

^F„:=B„y„=(B^-V^i^^zAy„, where A„=T^^, 
dt ~{ ^ ~ dt j ~[ z - 

and consequently systems of Schlesinger equations for the elements of Anj (or the coeffi- 
cient functions evaluated at Zj) are given in 118.6718.69^ . It is quite natural that systems 
governed by regular semi-classical weights preserve the monodromy data of the solutions 
Yn about each singularity Zj with respect to arbitrary deformations. 

Another theme we wish to develop is the evaluation of the above Toeplitz determinants 
in terms of generalised hypergeometric functions. Given a partition k = (ki, K2, • ■ • , K]v) 
such that Ki > K2 > ■ ■ ■ > i^N > one defines the generalised, multi-variable hypergeo- 
metric function through a series representation |46II81| 



(1.10) pF^^\ai, . . . ,ap;bi, . . . ,bq;ti, . . . ,tN) = ^ 



for p,q £ Z>o. Here the generalised Pochhammer symbols are 

AT 

(1.11) \4^'> ■■=X\^a-j + l)^^, 
the hook length is 

(1.12) n Hhj)+Khj)+^, 

where a{i,j),l{i,j) axe the arm and leg lengths of the (i,j)th box in the Young diagram 
of the partition k, and SK{ti, . . . ,tjv) is the Schur symmetric polynomial of A'^ variables. 
The superscript (1) distinguishes these functions from the single variable N = 1 functions 
and also indicates that they are a special case of a more general function parameterised 
by an arbitrary complex number d 1. 

In Section 2 we derive systems of differential-difference and functional relations for 
orthogonal polynomials and associated functions on the unit circle for a general class 
of weights and formulate the Riemann-Hilbert problem. In Section 3 we specialise to 
regular semi- classical weights and derive bilinear difference equations. In addition we 
calculate the deformation derivatives of the orthogonal polynomial system, derive a system 
of Schlesinger equations and show the deformations are of the isomonodromic type. The 
foregoing theory is utilised in the simplest case of three singularities and the A-recurrences 
derived in Section 4. In Section 5 the connection of the orthogonal polynomial theory 
with the Okamoto r-function theory is established. Application of our recurrences to the 
physical models described previously are considered in Section 6. 

2. Orthogonal Polynomials on the Unit Circle and Riemann-Hilbert Problem 

We consider a complex function for our formal weight w{z), analytic in the cut complex 
z-plane and which possesses a Fourier expansion 



(2.1) wiz) = Yl ^'^z'' »fc = / T^MOC 



-k 
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where T denotes the unit circle \(^\ = 1, with z = e'", 6* G (— t, tt]. Hereafter we wiU assume 
that z-'w{z),z-'w'{z) £ L{T) for aU j G Z. The doubly infinite sequence {wk}T=-oo a^re 
the trigonometric moments of the distribution w{e^^)d6 /2tx and define the trigonometric 
moment problem. Define the Toeplitz determinants 



:= det 



-MOC 



jj 2'rtiC, 
: det [w_,+j_fe](,<^. 



0<j,fc<n-l 



(2.2) =1/ n^-(OKr n ic.-ar, 

-^T" ,=1 ^'T*'.; i<j<fc<„ 

where e = 0, ±1. The last equality shows the determinants are equivalent to the CUE 
averages defined earlier and their generalisations. In certain circumstances the weight is 
real and positive, w(z) = 111(2) where the bar denotes complex conjugate, for example when 
all the monodromy parameters are real and certain independent deformation variables 
t £T, and so the Toeplitz matrix I^[w\ is Hermitian (see the example of the weight in 
section |1J but in general this will not be true. 

Furthermore consider the system of orthogonal polynomials {<^n(2:)}„gz>g defined with 
respect to the weight ■w{z) on the unit circle, assuming that none of the In[vj\ vanish. This 
system is taken to be orthonormal 

(2.3) / -^»(C)<^„.(C)04O=5-,n 

and the leading and trailing coefficients are defined by 

n 

(2.4) (t>n{z) = K,„Z" + InZ"^^ + TUnZ"^'^ + . . . + (t)n{0) = ^ Cn,j z\ 

j = 

where Kn is chosen to be real and positive without loss of generality. We also define the 
reciprocal polynomial by 

n 

(2.5) rn{z) := 2>41^) = E 

The orthogonal polynomials are defined up to an overall factor by the orthogonality with 
respect to the monomials 

(2.6) 0<j<n-l, 
whereas their reciprocal polynomials are similarly defined by 

(2.7) /'^„(^)^;(OCT=0 l<j<n. 

The system is alternatively defined by the sequence of ratios r„ = 4>n{Q)/K,n, known 
as refiection coefficients because of their role in the scattering theory formulation of OPS 
on the unit circle, together with a companion quantity fn (notwithstanding the notation, 
only when w{z) is real does fn equal the complex conjugate of r„). From the Szego theory 
|45| r„ and f„ are related to the above Toeplitz determinants by 

(2.8) r„ = (-l)":|M, f„ = (-l)"^"'t'"l 
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The Toeplitz determinants of central interest can then bo recovered from 

Further additional identities from the Szego theory that relate the leading coefficients back 
to the reflection coefficients are 

(2.10) 4 = 4-i + |</'n(0)|', 

, n— 1 

(2.11) ^ = Y,r,+,r,, 

(2.12) ^ = 

j = l 

Some useful relations for the leading coefficients of the product of a monomial and an 
orthogonal polynomial or its derivative are 

' ~ "I r Vn-l[Z) + %n-2 



Z 4>n{z) = 9n+2{z) + <Pn+l{z) 



Kn+2 \Kn+l lin+2 /tn+1 

(2.13) _|_ I ^"+1 ( _ _|_ mn _ mn+2 

\ Kn+1 \ Kn+2 Kn J «n Kn+2 



(j)n{z) + ■Rn-1 



(p'„{z) = n-!^4>„-i{z) + 7r„_2 

Kn—l 

Z<p'n(z) = nct>„{z) ^—(l)n-l{z) + 7r„-2 

f^n—l 

Z^(l>'n{z) = n (j)„ + liz) + I (n - 1)— - n ^""""^ I (l>„{z) + TTn-l 

where ' denotes the derivative with respect to z and where vTn denotes an arbitrary poly- 
nomial of the linear space of polynomials with degree at most n. 

Fundamental consequences of the orthogonality condition are the mixed linear recur- 
rence relations 

(2.14) K„<j)n+l(z) = K„+lZ(f>„(z) + (f>„+l{0)(l>l{z) 

(2.15) /tn0n+l(«) = K„+l(t>n{z) + <j>„+l{0)z4'„{z) 

as well as the three-term recurrences 

(2.16) Kn4>niO)(l}n+l{z) + Kn-l<i)n + l{0) Z<i>„-l{z) = [Wn (/>n+l (0) -|- + ! <^„ (0)2:] <^„ (z) 

(2.17) Kn(j)n{0)(j)n+liz) + «:„- 1 0„+l (O)«0* _ i («) = [K„(j>„+l{0) Z + Kn+1 0n (O)]0ji («) 

The analogue of the Christoffel-Darboux summation formula is 



(2.18) E'^i(^)'^^(C) = 



1 - < 



^2 ^ 'l>n+liz)<j)l+j^iO - (j}r^+l{z)(^>r^ + l(() 

1-zC 
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for zC + 1- 

Equations H2.16I2.17|I being second order linear difference equations admit other lin- 
early independent solutions i/)„(2:), and we define two such polynomial solutions, 
the polynomials of the second kind or associated polynomials 

(2.20) Vn(2) := /":5^^«;(C)[0n(C)-<^„(z)], n>l, := 1, 
and its reciprocal polynomial The integral formula for -i/;^ is 

(2.21) ^l{z) [ + ' wiOlz-MC ~ 4>Uz)], n>l, ^o*:=l. 

A central object in the theory is the Caratheodory function, or generating function of the 
Toeplitz elements 

(2.22) F{z) := / ^^w{0 



1^ 2'KiC, C, — z 

which has the expansions inside and outside the unit circle 



(2.23) F{z) = 



l + 2Er=i»fe^*. if|2|<l, 

-l-2E^=l^"-fe^"^ ifi^i>i. 



Having these definitions one requires two non-polynomial solutions en{z),en{z) to the 
recurrences and these are constructed as linear combinations of the polynomial solutions 
according to 



(2.24) e„{z) ■- ^„{z) + F(z)<j,4z) = [ 

J J 



dC, C + z 

■™(C)</'n(C) 



27rjC C,- z 



eliz) := V:(2) - F{z)cbUz) = -z" / ^ y wiOMO 

Jf Z-KiQ Q~ z 

Theorem 2.1 (^,[5^1,1151,201). il)n{z),iljl{z) satisfy the three-term recurrence relations 
12.1 (A and along with e„(2;), e* (2:) satisfy a variant of \2. 1412. 131 namely 

(2.26) «;„e„+i(z) = K„+ize„(2:) - (0)6^^(2:) 

(2.27) K„e*+i(2) = «;„+ie* (z) — <^„+i(0)ze„(z) 

Theorem 2.2 (|24,). The Casoratians of the polynomial solutions (j)„,(l>n,'<Pn,''pn o.'re 

(2.28) 0„+i(z)V'n(z) - i'n+i{z)(f>„{z) = (j}„+i{z)e„(z) - 6,1+1 (z)<^„(z) = 2- 



/>n + l(0). 



4n-Fl(0) „+l 



(2.29) </>;^ + l(z)V'n(z) - V'n + l(2:)0n(z) = <i>*n + \{zYn.{A - 4 + l(2)C(2:) = 2- 

(2.30) (l)„{z)ip*{z) + jp„{z)(f>*^{z) = (^„(z)e* (2) -f- e„(2)(^* (2) = 22" 

We will require the leading order terms in expansions of <j)n{z), (t>^{z), £^(2), e^(z) both 
inside and outside the unit circle. 
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+ 7 ""^^^"^ ^ ""^^^"^ "^^ ) z-" + 0(g-3) \z\>l 



0n + l(O) ^„ + l ^ / nl, <^n + 2(0) (^n + l(0) + 1 \ „+2 | Q(-;'^ + 3-) |2|<1 



Corollary 2.1. T/ie orthogonal polynomials (j)„{z), (p^{z) have the following expansions 

( MO) + -^{lin<Pn~l{0) + MO)L-l)z + 0{Z^) \z\ < 1 

(2.31) Mz) = < ''"-1 

fK„ +r„2 + 0(2:^) |z| < 1 

(2.32) = < - 1 - 

\MO)z" + {K„Mi{0) + MO)ln-i)z"-' +0{z"-^) \z\>l 

whilst the associated functions have the expansions 
(2.33) 

2 ^"^'^''-*) <?!.„+i(0)^_i , ( kI <?!>„+2(0) 9i„+i(0) 
(2.34) 

^e„(2;) 

Z ■ 1 _ _^ j -"^-'"^^ _ iriLLl ) + o(z-3) |z| > 1 

The z-derivatives or spectral derivatives of the orthogonal polynomials in general are 
related to two consecutive polynomials 27 and we generalise this with the following 
parameterisation. 

Proposition 2.1. The derivatives of the orthogonal polynomials and associated functions 
are expressible as linear combinations in a related way (' :— d/dz), 

(2.35) W(z)Mz) = Qr^{z)Ml{z) - {^u{z) + V (z))^^) 

(2.36) W{Z)^*:{Z) = -e;(2)<^: + i(2) + {K(Z) - V(z))rn{z) 

(2.37) W{z)e'„{z) = e„(2)e„+i(2) - (^^(z) - V{z))e4z) 

(2.38) W{z)e*:(z) = -e;(z)e;+i(z) + (Kiz) + Viz))eUz) 
with coefficient functions W {z),V (z) independent of n. 

Proof. The first, 12.351 . was found in 1271 where the coefficients were taken to be (their 
notation An, B„ should not be confused with our use of it subsequently) 

, . . _ /s:n-10n+l(O) 2:9n(z) 

^ ' ^„0n(O) W{Z) 

(2.40) Bn ^ ^ («„(.) + Viz) - + e„(.)) . 

The other difi^erential relations can be found in an analogous manner. □ 

The coefficient functions &ri{z),Q,n{z),Qn{z),^n{z) satisfy coupled linear recurrence 
relations themselves, one of which was reported in |27| . The full set are given in the 
following proposition. 

Proposition 2.2. The coefficient functions satisfy the coupled linear recurrence relations 

(2.41) n„(z) + nn-,{z) - ( + ^z\ Qn{z) + {n- 1)^ = 
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, K„</)n+2(0) , , K„^icj>n+l{0) 0n+l(O) W(z) _ ^ 



(2.43) + n„_i(2:) h -^-r—z Q^{z) - n = 

\ K,i (pn[U) J Z 



(t>n[0) J 

Kn + l<Pn + l(0) Kn(pn(0) Kn Z 



(2.45) n„+i(2) + n:(^) - f + + ^(^e„(^) - e:(^)) = 

(2.46) n„iz) - n„+,iz) + ^(z + iji±M^^ji±m) e„+,(.) 

+ Mi±iM0;(,) _ !^.e4.) - ^ = 

(2.47) ^UAz) + f^„(z) - f ^ + 1^^) e:+i(2) 

(2:e„(2;) - Q„{z)) = 

(2.48) n:(.) - n:+i(^) + ^ f 1 + ^niiM^iiiiM,) e:^,(,) 

H z<dn\z) B„(2) = 

Proof. The first (ITITI was found in |27| by a direct evaluation of the left-hand side using 
integral definitions of the coefficient functions, however all of the relations follow from the 
compatibility of the differential relations and the recurrence relations. Thus 12.4112.421 
follow from the compatibility of llO^ and lITTSl . ^T^ . ^TMi follow from ^TM and 
(ITTTIl . il4Sll>.46ll follow from the combination of llia5ll>.:j6ll and lITTIIl . and ^lA'lMA^X 
follow from the combination of 12.3512.3611 and 12.151 . □ 



Remark 2.1. The relations given above are obviously not all independent, as for example 
we note that 12. 4H can derived from 12.451 1 with the use of 12.5(JI I below. 



11 



Corollary 2.2. Some additional identities satisfied by the coefficient functions are the 
following 



(2.49) 

(2.50) 
(2.51) 



l(0) 



0n(O) 



Quiz) 



zQn^l{z) = - Z@^(Z) 6n-l(2) 



MO) 



Q.l{z) - n„(2) = -^!ii±i(ze„(z) - e:(z)) + n 



W{z) 



n*„{z) + n4z) = 1- 



<^n + l(0)</>„ + l(0) 



(^n + 2(0) , , Kn+1„*/ n 
6n + l(2) H 0n(^) 



(0) 



+ 



W{z) 



For a general system of orthogonal polynomials on the unit circle the coupled recur- 
rence relations and spectral differential relations can be reformulated in terms of first order 
2x2 matrix equations (or alternatively as second order scalar equations). Here we define 
our matrix variables and derive such matrix relations, and this serves as an introduction 
to a characterisation of the general orthogonal polynomial system on the unit circle as the 
solution to a 2 X 2 matrix Riemann-Hilbert problem. 

Firstly we note that the recurrence relations for the associated functions e„(2:),e*(z) 
given in 112. 2612. 27|l differ from those of the polynomial systems 112.1412.15^ by a reversal 
of the signs of (j>niO), <^n(0). We can compensate for this by constructing the 2x2 matrix 



(2.52) 



Y„{z) 



4>n{z) 



w(z) 
en{z) 



w{z) I 



and note from II2.29I I that dety„ = ~2z^ lw(z). 



Corollary 2.3. The recurrence relations for a general system of orthogonal polynomials 
i2.i4t2.rEli) and their associated functions t2.26\2.FJ\) are equivalent to the matrix recur- 
rence 



(2.53) 



Yn + l :— M„Yn — — 



1 / K„ + lZ (t)n + l{0) 

lin \(j>„+l{0)z K„+l 



ith according to 12.1(A) , det M„ — z. 



Corollary 2.4. The system of spectral derivatives for a general system of orthogonal 
polynomials and associated functions 12.3512. 3lA) are equivalent to the matrix differential 
equation 



Y' — 4 V 



(2.54) = 



Wiz) 



n„(z) + v{z)~'^ze„{z) 



l(0) 



Pn+l 



(0) 



V 



zOn{z) 



il*(z) 



V{z) 



Quiz) 



Proof. This follows from llij.;j5l!^.38t and employing lli^.l4l:i.l5l!^.i!6li>.ij7ll . 
Remark 2.2. Compatibility of the relations 112.5311 and 12.541 leads to 
(2.55) Af; = A„+iM„-M„yl„, 



e;(z) 



□ 
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and upon examining the 11-component of this we recover the linear recurrence I2.46|l . the 
12-component yields 1)2.45^ . whilst the 21-component gives l|2.47|l and the 22-component 
implies 12.481 . 

Remark 2.3. There are, in a second-order difference equation such as l|2.26|l or 12.2711 . 
other forms of the matrix variables and equations and these alternative forms will appear 
in our subsequent work. Defining 



(2.56) X„(z;t):= 



/ , / N gn+l(z) ' 
/ 4>,i + l{z) -— 



(t>n{z) 



X*{z;t) ■- 



■w{z) 



Vn+l{Z) 



we find the spectral derivatives to be 



(2.57) W{z)X'^ = 

V e„(z) 



7Tn-2e„+i(z) \ 

«:n + l(Pn + l(0) Xn 

-^^{z)-V{z) I 



(2.58) W(z)X: = 



-J7„(2) - \/(z) + (n+ 1)— — y 77^^'3"+i(^) 1 V* 

2 K„ + 10„ + l(O) I A„. 



^e:(2) 



Another system is based upon the definition 
(2.59) Z^{z;t) 



I (?!>„+i(z) -— 

w{z) 

<Pn{z) 



w{z) 



and in this case the spectral derivatives are 
(2.60) W{z)Z'„ 

(-ni{z) - v{z) + ^^e;(2) + (n + 1) 



<^n + l(0) 



e:(2) 



r!;(2)-i^(2) 



<Pn + Az) -f-^ > 

w;(2) 

, . s en(2) 
. (I>n{z) ——- I 

\ -wiz) / 



«n<^n + 2(0) 
^n+1 



e„+i(2) 



n;(2)-v'(2)-^e;(2) 



(2.61) 

w{z)z: 



n„{z)~V{z)^ ^^zQ„{z) 



i(0) 



V 



Qn{z) 



z'^Q^+iiz) 



^n+1 



-n„{z)-v{z) + 



-29,1(2) 



Kn+l ' ' ' ' ' ' Hn+l ' I 

We end this section with a characterisation of a general system of orthogonal poly- 
nomials on the unit circle (and their associated functions) as a solution to a particular 
Riemann-Hilbert problem. 

Proposition 2.3. Consider the following Riemann-Hilbert problem for a 2 x 2 matrix 
function Y : C ^ 51/(2, C) defined in the following statements 

(1) Y{z) is analytic m {2 : \z\ > 1} U {2 : \z\ < 1}, 
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(2) on z G E where E is the oriented unit circle in a counter- clockwise sense and 
+ (— ) denote the left(right)-hand side or interior(exterior) 



1 w[z)/z 
1 



(2.62) Y+{z) = Y^{z) 
(3) as z ^ oo 

(2.63) Y{z)={l + 0[z--))('^^^^ 



(4) as z 

(2.64) Y{z) = {l + 0{z)) 



/0(1) 0{z-') 
lo(l) O(z-) 



It is assumed that the weight function w{z) satisfies the restrictions given at the beginning 
of this section. Then the unique solution to this Riemann-Hilbert problem is 

I <fi„{z) enjz) \ 



(2.65) Y{z) 



^) 



n> 1. 



Proof. We firstly note from the jump condition 12.6211 that Vn, Y21 are entire z G C. From 
the 11-entry of the asymptotic condition 12.61311 it is clear that Y\\ — ■Kn{z) a polynomial 
of degree at most n. Similarly Y21 = o-n{z) from an observation of the 21-component. 
From the 12- and 22-components of the jump condition we deduce 

(2.66) Y+^2~Y.^2 = ^Y^i, Y+22-Y-22 = ^Y2i, 

z z 

and therefore 

(2.67) Y,2 = / J^!^!(CWC) , y,, = ^ J^MOMC) 



27ri(^ ~ z Jy 2iiiC, C, — z 

Consider the large z expansion of Y\2 implied by the first of these formulae 

(2.68) Y,2 = -z-^ ^ _^^(0^„(C) + 0{z~^). 

According to 12.6311 the integral vanishes and so ir-niC,) is orthogonal to the monomial 
Now take the small z expansion 

(2.69) F12 = / /^^(C)^n(C)C^+0(z"-^). 

From the 12-component of the condition 12.6411 we observe that all terms in the sum vanish 
and we conclude the 7r„((^) is orthogonal to the monomials C,, . . . , and the first term 
which survives has the monomial C". Thus Tin{z) oc (j)n{z), and from the explicit coefficient 
in the 11-entry of 12.6311 iin{z) is the monic orthogonal polynomial <j}„{z)/K„. We turn 
our attention to 122 and examine the small z expansion 

(2.70) Y22 = E / A^^iCWiOC^+oiz"). 

The 22-component of Ij2.64|l tells us that all terms in the sum vanish and consequently 
o"n(C) is orthogonal to all monomials ...,<^". Therefore o-„{() oc (j>^{z) and we can 



14 



determine the proportionality constant from the 22-component of the asymptotic formula 
and comparing it with 

(2.71) Y22 = ~z~^ f T^wiOar^iO + 0{z-''), 



If 27ri(^ 

to conclude Un{() = «„<;!)* (z). Finally we note that 

f rfC HOMO _ 1 .^ f dC wiQK{^__l 

^ ' AaTTiC C-z '2z"^^' Jj2n<: C-z ' 2z"^'' 

when n > 0. We also point out det Y — —z"^^. □ 

Remark 2.4. Our original matrix solution ¥„ specified by 12.521 is related to the solution 
of the above Riemann-Hilbert problem by 

1 




K„/ \ w{z)/ \ u K„y y 22 

Our formulation of the Riemann-Hilbert problem differs from those given in studies con- 
cerning orthogonal polynomial systems on the unit circle with more specialised weights 
m. We have chosen this formulation as it is closest to that occurring for or- 
thogonal polynomial systems of the line JISJ, the jump matrix is independent of the index 
n which only appears in the asymptotic condition and it is simply related to our matrix 
formulation H2.52|l . 

3. Regular Semi-classical Weights and Isomonodromic Deformations 

All of the above results apply for a general class of weights on the unit circle but now 
we want to consider an additional restriction, namely the special structure of regular or 
generic semi-classical weights. 

Definition 3.1 (gZI)- The log-derivative of a regular or generic semi-classical weight 
function is rational in z with 

(3.1) W{z)w'{z) = 2V{z)w{z) 

where V{z),W{z) are polynomials with the following properties 

(1) deg(W) > 2, 

(2) deg{V) < desiW), 

(3) the m zeros of W{z), {zi, Z2, ■ ■ ■ , Zm} are distinct, 

(4) the residues pk = 2V{zk)/W'{zk) ^ 1. 

The terminology regular refers to the connection of this definition with systems of 
linear second order differential equations in the complex plane which possess only isolated 
regular singularities, and we will see the appearance of these later. An explicit example 
of such a weight is that of the form 111(2) = W^^^(z — ZjY^ with Zj 7^ z^ for j 7^ fc, which 
are also known as generalised Jacobi weights. In addition we will assume the polynomials 
defined above take the following forms 
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The above definition is restrictive and has been generahsed by relaxing some of the con- 
ditions in a series of works |41) . |89| . |4()| . In these works the orthogonal polynomial 
systems were characterised by integral representations of semi-classical linear functionals 
with respect to certain paths in the complex plane. 

It follows from these definitions that the Caratheodory function satisfies an inhomo- 
geneous form of Ijl-i.l^ . 

Lemma 3.1 (^lEZl)- The Caratheodory function 12. satisfies the first order linear 
ordinary differential equation 

(3.3) W{z)F'{z) = 2V{z)F{z) + U{z), 

where U{z) is a polynomial in z. 

This lemma leads to the following important result. 

Proposition 3.1. The coefficient functions 0„(2:), 0* (z), f2„(2;), Jl* (z) are polynomials 
in z of degree m — 2,m — 2,m — l,m — 1 respectively. Specifically these have leading and 
expansions of the form 



(3.4) e„(^)^in + l + Y^p,)-^z"^-' 



J = l 



{iiL III, iiL lit, ^ I / r\\ 

-[{u+I + Y: p.) E - E ^^^^1^ + (n + 2 + ^ p,) , ° 



(3.5) e„{z)^[2V{0)-nW'{0)]- 



i^n + T p^.)^!l±iMlM _ 2^1^ \z^-' + 0{z 
.(0) 



m — 4\ 



/>n + l(0) 

+ {[2V'{0) - ^nW"m-^^ + [21/(0) - (n - l)vy'(O)] ^^^^^^ 

1^ (f>n + l(0) K„-l(pn + l{0) 

+ ([{n + l)VK'(O) - 2Vm^ - [(n - l)Vy'(O) - 21^(0)]-^) + 0(^' 

(3.6) e:{z)^-in + f:p,)J^.-^ 



+ + E p^) E - E ^^-^^ix^ + (- + 1 + E 



j = l + ! (0) ^ (^n + l(0) «:„ + ! 

^ 1 , ^ '^n(^n-l(0) + <^n(0)/n-l 1 m.-3 
-(71-l+2^Pj) 



■^" + 1(0) 



(3.7) Ql{z) = -[2^(0) - (n -f l)IV'(O)]^ 



+ I - [21/'(0) - i(n + 1)VK"(0)]-^ - [21/(0) - nW/'(0)]-^ 

+ \{n + 2)T^'(0) - 2\/(0)] f - Jlj^h^±l\ L + O(z^) 
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(3.8) J7„(2) = (1 + V2E^^)^'""' 



3 = 1 

{mm m m 

3 = 1 3 = 1 3 = 1 3 = 1 

K,i+2/tn+l 0n + l(Uj Kji + l J 



(3.9) J7„(2) = y(0) -7iiy'(0) 

+ {v'm - inVy"(0) + fl/(0)^ + [V-(O) - nw'm]^) 



+ [1/(0) -nW'(O)]— - [V(0) - (n+ l)VK'(0)]-^L + O(2'') 



(3.10) n:(2) = ~i/2^p,z™-i 



3 = 1 
m m 



I ^1 ^1 ^1 ^/''Kn + l<^„ + l(0) + 

+ 0(2"-^) 

(3.11) n:(2) = (n+l)W^'(0)-l/(0) 

+ 1 1 (n + 1)W^"(0) - l/'(0) + [(n + 2)W'(0) - 2V(0)] ^iiilM _ v[/'(0)i^ ] 

L «;n+2Kn+l (Pn+l(0) Kn+1 J 

+ 0(.^) 

Proof. Following the approach of Laguerre we write F{z) in terms of (l>n{z),ipn{z), en{z) 
and use (13.311 to deduce 

(3.12) Q = WF' ~2VF -U 

(3.13) =^f '"r'^" V"2V^ '"T'^" ~U 



z 



2 



(3.14) ^ " v-nv-n; , ^ v v^n l^V^i ^ l[y _ 



The numerator of the first term is independent of e„, and so is a polynomial in z, and we 
denote this by 

(3.15) 2^Ii±iMz"e„(2) = W(-Ke'„ + e^^'^) + 2T/0„e„. 

Given that this is a polynomial we can determine its degree and minimum power of z by 
utilising the expansions of e„ both inside and outside the unit circle, namely 12.3112.331 . 
We find the degree of the right-hand side is n + m — 2 so that Qn{z) is a polynomial of 
degree m — 2. Developing the expansions further we arrive at l|3.4|l . An identical argument 
applies to the other combination 

(3.16) 2^^i±i^2"+^e;(2) = W{^le': - <<^:') - 2V4>lel 
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and On{^) is also a polynomial of degree m — 2 with the expansion 

lEHJ- To establish 

H,'-{.8|l we utilise the other form of 0^(2:) and l|2.28|l to deduce 

(3.17) W{iP„(l}'„ - (/.nO + 2V0„7/)„ - = 2^^^±^z"e4z) 

(3.18) = [(t>n + l1pn - 'll)n+l4>n]Qn{z). 

Separating those terms with (j)n and as factors we have 

(3.19) {e„(z)<?i„ + l - W(t)'r. - V(l>n] = {e„(2)^/>„ + l - W i^'^ + " f^'/'n } (/-n , 

so that this polynomial contains both 0„ and ■i/'n as factors and can be written as ^n4>ni>n 
with Q.n{z) a polynomial of bounded degree. This latter polynomial can be defined as 

(3.20) 2 '^" + ^(°) z"»„(z) = W{i^^ + l<f}'„ - 0n + l^/';) + + l + Vn-^n + l) ^ U <i>n<l>n + l 

(3.21) = W{€„ + l(p'„ — (pn + ie'^) + V{(l)n<:n + 1 + Enl^n + l). 

Again employing the expansions 12.3112.3311 we determine the degree of Qn (z) to be m — 1 
and the expansion 13.81 follows. Starting with the alternative definition of Q^{z) and 

(3.22) vK(0:c - c<^;') - 2vr„rr. - u<t>:^ = 2^^^±^z"+'e:{z) 

(3.23) = - Vn + 10n]©n(2)- 

and using the above argument we identify for the polynomial Q^{z) 
(3.24) 

20riJ-l_(O) ^n + lj^. j.^^ = W(-?/)* + l(?!)*' + (I>n + l1pn) ~ V {(f)*„-tp'^ + i + V'n<^n + l) " f7(?!)* (/)* +i 

(3.25) = V(^(-e^+i<;/)*i' + (;/)*+ie*') - V'(<;/)*e*+i + e^fji^^^). 

The degree of f2* (z) to be m — 1 and has the expansion 13.1UII . □ 

Remark 3.1. Solving for (f>'^ and between 13.1511 and 13.211 leads to 12.351 and 12.371 . 
whilst solving for </!>*' and using 13.1613.251 yields 12.361 and 12.381 . 

Furthermore, in the case of a regular semi-classical weight function, the matrix A„{z; t) 
has the partial fraction decomposition 

A„j (t) 



(3.26) A„{z;t)~Yl 



under the assumptions following 13.1|l and the residue matrices are given by 
(3.27) 

A - Pi 



2V[z 



-Q,n{zj) — V{zj) H ZjQn{zj) — Qn{zj 



.... riX^jl 



\ ZjQ„(Zj) n„{Zj)~V{Zj) Qn{Zj)j 



Using the identity 12.501 we note that TrA„j — —pj and TrA„(2;;t) — —w' (z) /w{z). 

Bilinear residue formulae relating products of a polynomial and an associated function 
evaluated at a singular point will arise in the theory of the deformation derivatives later 
and we give a complete list of results for these. 
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Corollary 3.1. Bilinear residues are related to the coefficient function residues in the 
following equations 



(3.28) 
(3.29) 
(3.30) 
(3.31) 
(3.32) 
(3.33) 
(3.34) 
(3.35) 
(3.36) 
(3.37) 



Kn Z V yZj ) 



<t>niZj)en + l{Zj) 
'^n+l(2j)en(2j) 



2V{z,) 

o 0n + l(O) . ) + V(2j) 

«;„ ' 2V{z,) 

n + 1 {Zj 



JL 

V{z, 
zl 



2V(z,) 

0.rXZ]) - V{Zj) - 



jQn{zj) 



V(Z, 



V{z, 



nl{z,)-V{z,)-'^Ql{z,) 

nn{Zj) + V{Zj) - '^^Zje„{Zj) 



ni(z^) + v{z,)^^e:{z,) 



Proof. These are all found by evaluating one of H3.15^ . I|3.16^ . I)3.21|l . or Ij3.25|l a.t z — Zj 
and using (j^HEUli- □ 
Remark 3.2. The initial members of the sequences of coefBcient functions {O„}^oi 
{e;}SLo, {^n}^=o, mn=o are given by 
(3.38) 

2V{z) - U{z) 



eo{z 

(3.39) 

ei(2 

(3.40) 
(3.41) 

el(z 

(3.42) 

Qo{z 
(3.43) 

Qo{z 



^kIz^{2V{z) - U{z)) - 2Kiflii(0).z(/(2) - 2Ki(?ii(0)W(2) - (j)i{Q){2V{z) + U{z)) 



2V(z)-U(z) 



^<j}l(Q)z^{2V{z) - U{z)) - 2ki<^i(0)2[/(z) - 2ti-,(j}i{Q)W (z) - «:?(21/(z) + U{z)) 



z[2V{z)~U{z))~^ 



201 (0) 



- -{2V[z) + Uiz))^''^'^ 



2(l>i{Q)z' ' ' ^ " 2 
One can take combinations of the above functional-difference equations and construct 
exact differences when z is evaluated at the singular points of the weight, i.e. W(z) = 0. 
The integration of the system is given in the following proposition. 
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Proposition 3.2. At all the singular points Zj,j — 1, . . . m, with the exception of Zj — 0, 
the coefficient functions satisfy the bilinear identities 



(3.44) 
(3.45) 

(3.46) 
(3.47) 



«:n + l<?!>n + l(0) 



^.2 ^„</>„+2(0) 



2 , i2 



J7„_i(z,)-^%i#e4^.) 



2je„(2j)e„+i(2j) + V (zj) 
zie;(zj)e;+i(2j) + v'^izj) 



kI 0n(O) 



e„(2,)e:_i(2,) + \/^(2,) 



(3.48) 

(3.49) 



«n-l <?!'n + l(0) 



1^1 0n(O) 
0n + l(O)<^„ + l(O) 



2^9; (2j) 



Kn-l'/'n + l(0)<^n(0) 2p,. / 

2jO„^2jje'n 



-1 U:.) 



en(2,)e:(2,) + \/2(2,) = 



(2j) - ■^^2je„(2j) 



n*„{zj) 



Hn + l 



e;(2^) 



First Proof. We take the first pair of identities H3.44^ and H3.45t as an example for our 
first proof. Multiplying the Q„, Q„-i terms of II2.41|I by the corresponding terms of H2.42^ . 
evaluated at a singular point z = Zj, one has an exact difference 
(3.50) 



r)2 f \ «ra<^n + 2(0) ri ^ \n f \ ^n-l<^n + l(0) „ f f \ 

S2„„i(2j) — 7 7-:;ZjC>n(Zj)<3n + l{Zj) , . Zjt)n^l(Zj)<dn(Zj), 



^^y,^..,.,,^..^.y.j, /.„<^„(0) 

assuming none of the Zj coincide with —r„+i/r„ for any n. Upon summing this relation 
the summation constant is calculated to be 

«:o</>2(0) 



(3.51) 



0.l{zj) 



zjQo(zj)ei{zj) = V (zj), 



by using the initial members of the coefficient function sequences in 13.4213.3813.391 . The 
result is (13.441 . whilst the second relation follows from an identical argument applied to 
j2.43l2.44t . □ 



Second Proof The three pairs of formulae 13.4413.451 . 13.4613.471 and 13.4813.491 arise 
from the fact that at a singular point Zj the determinant of the matrix spectral derivative 
must vanish. Thus 13.441 and 13.451 express the condition that the determinant of the 
matrix on the right-hand sides of 12.571 and 12.581 vanish respectively. It can be shown 
that the same condition applied to the right-hand sides of 12.611 and 12.601 implies 13.461 
and 13.471 respectively when one takes into account the identities 12.461 . 12.491 . 12.501 
and The last pair are a consequence of det{W A„(zj;t)) = along with the identity 

Third Proof. All the bilinear identities in Proposition 13.21 can be easily derived from the 
residue formulae 13.2813.371 by multiplying any two of the above formulae and then fac- 
toring the resulting product in a different way. Thus 13.441 arises from multiplying 13.301 
and 13.311 and then factoring the product in order to employ 13.281 . Equation 13.461 
comes from multiplying 13.281 and 13.291 with n n — 1, using the recurrences 12.151 . 
12.271 with n n — 1 to solve for (j)^_i{zj),e^_i{zj) and employing 13.281 along with 
and Km setting n 71 — 1. Equation 13.481 is derived by multiplying 13.281 
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and HH.29|I and then factoring using lj8.H4|l and The reciprocal versions follow from 

similar reasoning. □ 

Remark 3.3. It is clear from the first proof that the bilinear identities given in Proposi- 
tion 021 can be straightforwardly generalised to ones that are functions of z rather than 
evaluated at special z values. They can be derived directly from Proposition 12.21 so ap- 
ply in situations where the weights are not semi-classical, and contain additional terms 
with a factor of W{z) and sums of products of other coefficients ranging from j — 1, . . .n. 
However because we will have no use for such relations we refrain from writing these down. 

Remark 3.4. If 2; = is a singular point then the limit as z ^ may be taken in the 
product of 12.4112.4211 ■ however this does not lead to any new independent relation but 
simply recovers 

nn(0) = ViO)-nW'{0). 

We now consider the dynamics of deforming the semi-classical weight 13.21 through a 
i-dependence of the singular points Zj{t), 

m 

(3.52) 

where ' := djdt. Given this motion of the singularities we consider the t-derivatives of the 
orthogonal polynomial system. 

Proposition 3.3. The deformation derivative of a semi- classical orthogonal polynomial 



(3.53) = { ~ ^^Y.p^- + y^Y.p^-^r^-i^j)^u^j)^—}M^) 

whilst that of a reciprocal polynomial is 

(3.54) 4>l{z) = {- — + ^/2Y^p,^-^z]--el{z,)<t>r.{z,)^-]<t>l{z) 

m 

- {^'^^ P3—zj"aZ])<P*JZ])—^ ]<Pn{z) 

The deformation derivative of an associated function is 

m 

(3.55) Uz) = { - - - y2Y,pAzr^Uz,)Mz,)^—}er.{z) 

V hxi 77, Z'j Z Z j } 

I V2 y^Pj— Z^'"tn {Zj )(l>nizi) I < (Z) 

L Zi z ~ Zi > 



J = l 



+ 

J = l 
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and that of a reciprocal associated function is 

(3.56) kl[z) = I - — - ^/2^ pj^z]'"" en{zj)(l>l[zj)——\el{z) 

m 



+ 



Proof. Differentiating the orthonormality condition 

and using l)3.52|l we find 

kn f dC, , , • ■r-^ . [ dC, , . 1 



= —Si,o + / 7r~7™(^)'^"'?^"-^ ~ Yl Pj^i / ^r~7'^(^^)7 <i}n(t>n-i, i = 0,.. 

Kn J 2-KiC ^ J 2mC C ^ 

Now 

J 2^"'(^)c^'^"(^)'^"-'(^)=y 2^"'('^)'^"(^) c^z 



+ 1 



so that 



^--^..o+ '^"-;^^''^ 6„(z), n>0 
z 2z 



\ l^ri . Zj J ^.^^ Zj J ^TTJC, 

In addition we can represent (j>„-i{z) as 

n 

4>n-i{z) = E 5i^j(j)n-^j{z) 
3 = 

= i 2^^^^)'^— l^z ■ 

Writing the Kronecker delta in a similar way the whole expression becomes 

for all < j < n and 13.531 1 then follows. The second relation follows by an identical 
argument applied to 

^U.(C)0n-4C)0^=<5.,O. 
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The derivatives of the associated functions H8.55|l . foUow from differentiating the 

definitions 12.24|l . H2.25|l and employing the first two resuhs of the proposition along with 
the relation (lOHt . □ 

Corollary 3.2. The t-derivatives of the reflection coefficients are 

(3.58) fe-V.E/.|^^^-#^^ 
Proof. An alternative form to 1^3 .53^ is 

(3.59) <^„(c) = -{ — + Y.p^-] -^"(0 + y^Y^pA^niz,) <^„_,(^-^)<^„_kc), 

V^" J = l j = l 

and by examining the coefficients of (^^^ , we deduce that 



^--V^^l^P.-^. .n(^.)0„-r(..). 
i=i 

Noting that the derivative term of 13.2f H vanishes when z = zj and employing 12.2811 we 
arrive at 13.5711 . The second equation, 13.5811 . follows by identical reasoning. □ 

Sums of the bilinear residues over the singular points are related to deformation deriva- 
tives in the following way, 



^7^ = ~Ep^v+ y^Y.pA^r-r.{z,)<i>i{z,\ 

m 

(3.60) = -y2Y,p3-zr^U^j)M^3), 

X. ■ ■ m . 

<Pn(0) K,, 2j 1 / ^" / \ . / \ 

^ + - + Ep.- - E p.-.. enizM^,), 

fn r.^\ 0'l + l(O) P3 . _1„ , . 

xtS + - = -v^xTk E P.^^7"^:(%-)'^n(%) 

<^n(0) «:„ (^„(0) ^ 2j 

For the regular semi-classical weights we can also formulate the system of deforma- 
tion derivatives as a 2 x 2 matrix differential equation and demonstrate that the system 
preserves the monodromy data with respect to the motion of the singularities Zj{t). 

Corollary 3.3. The deformation derivatives for a system of regular semi-classical or- 
thogonal polynomials and associated functions 13.5513.561) are equivalent to the matrix 
differential equation 

(3.63) Yn~ BnYn = \Bo^^Y-^Anj\Yr,. 
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where 



(3.64) 



Bo 







(0)+K„<?i(0) 



V 



Proof. This follows from a partial fraction decomposition of the system 113. 5313. 5611 and 
using j8.6()l3.62t . □ 

In the case of the pair l|2.5H|l . Hi-i.63|l compatibility implies the relation 

(3.65) Mn = Bn + lMn " M„B„, 

however there are no new identities arising from this condition. Taking the 11-component 
of both sides of this equation we see that it is identically satisfied through the use of 
l|2.46|l and H3.62|l . Or if we take the 12-components then they are equal when use of made 
of and H3.6()I3.57|I . In a similar way we find both sides of the 21-components are 

identical when we employ (12.431 and 13.621 . Finally the 22-components on both sides are 
the same after taking into account 12.4811 and 13. 6013. 621 . 

For the pair of linear differential relations 12.541 . 13.631 compatibility leads us to the 
Schlesinger equations 



(3.66) 



^^^^ 



Again there is not anything essentially new here, that couldn't be derived from the system 
of deformation derivatives 13.5313.561 . but it is an efficient way to compute the deformation 
derivatives of bilinear products. Employing the explicit representations of our matrices 
Ak we find the following independent derivatives in component form 



(3.67) 

d pj 



dt 2V{zj 



Pj l<^»+i(0)|= 



■E- 



ZjQn{Zj J 



Zk Pk 
Zk 2V{zk) 



Pj <fi„+i{0) d 



2V{z, 



dt 



('^n<^7i(0))6„(2j 



[Zk0nizk)<d„(zj) - Z.jQ„{zk)'dnizj)] , 



(3.68) 



d Pj 0n + l(O) ^ 
dt 2V{Zj) Kn 



Pj (t)„+i{0) f k„ \ - ij - Z k Pk 



V{z, 



Zk 2V{zk) 



Qn{zk)[0.n{zj) ^^2j0„(2j)] — 0,i (zj ) [Sl„ (zfc ) ^^Zfc6n(2fc)] > 



Pj 4>n + l{Q) { kn 



V{Zj 



'-Zi&l(Zi 



+ 



d 



(K„(^„(o))[ri* (zj) 



Zk pk 



l^n(t>u + l{0) dt K„ f-'_ Zj - Zk2V{zk) 

2:fc8* (2:fc)[f2* (zj) - -^^^e* (zj)] - 2je*(2j)[f7* (zfe) - -^^^e* (zfe)] > 
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Remark 3.5. The fact that the deformation equations satisfy the Schlesinger system of 
partial differential equations should be of no great surprise as the isomonodromic proper- 
ties of the regular semi-classical weights are quite transparent. In the neighbourhood of 
any isolated singularity \z — Zj\ < 5 the Caratheodory function can be decomposed 

(3.70) F(z)=f,{z)+CM^), 

where fj (z) is the unique, holomorphic function in this neighbourhood and Cj a coefficient. 
The analytic continuation of ¥„ around a closed loop enclosing the singularity is easily 
found and defines the monodromy matrix Alj, 



(3.71) = V-„L^+,M„ 



1 Cj{l- e-^""'"^ 
e-^""'"^ 



From the definition 12.221 one can confirm that Cj = l/isin(7rpj), so that the Mj is quite 
naturally independent of the deformation variables Zj or t {pj being constant). 



4. The simplest Semi-classical Class: Pvi System 

Here we consider the application of the general theory above to the simplest instance 
of the semi-classical weight, namely m = 3 singular points with two fixed at z = 0,-1 
and the third a variable at z = — 1/t. 

Explicitly we consider the unitary group average (11.311 where t = e"*, L0i,uJ2,fJ. £ C 
{iJjLj — uJi i: iu!2) and ^ € C. In the initial formulation t £ T but will be analytically 
continued off the unit circle. The weight function 

(4.1) w{z)=t'''z-'"''{l + zf'^^l+tzf'' I ' , 

is known as a generalised Jacobi weight, with branch points at z = 0, —1, — 1/t, oo. When 
t £ T, n,u!i,u>2, £, € R and C < 1 this weight is real and positive. The Toeplitz matrix 
is then hermitian and as a consequence f„ is the complex conjugate of r„, but generally 
this is not the case. The Toeplitz matrix elements can be evaluated in terms of the Gauss 
hypergeometric function, and there are several forms this can take which exhibit manifest 
analyticity at either of the special points t = 0,1, oo. 

Lemma 4.1. The Toeplitz matrix element Wn for the weight {4- Ij , under the restriction 
K(/i), 5R(a;i) > — V^; given in terms of hypergeometric functions analytic at t = 0,1 

(4.2) t^Wn = p., ^ ^ T'^^VT/^ —r2Fi{-2ti, -n-tJ.-Lo;l-n-ti + ij;t) 

1 (1 + n + fj, + Lojl (1 — n — fi + ijj) 

£. ±7ri(n + M-c:') r(2^ + l)r(2LJl + 1) 

_l_n--\-fj. — uj / ^\2fj.-\-2uJi-\-l 

2^1 r(2p + 2lui + 2) ^ ~ ' 

X 2-Fi (2^ + 1, 1 + 71 + ^ + oj; 2/i + 2u;i + 2; 1 - f) 
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where the ± sign is taken accordingly as Im(t) ^ 0. This can also be written as 



(4.3) fwr. = n + e 



■i 



p±7ri(ri + /i — Lj) 



r(2a;i + 1) 



2i sin7r(n + /i — lj) J r(l + n + + (jj)r(l — ti — + lj) 

X 2-Fi( — 2/1, —71 — /i — — n — /i + Lj;t) 
^±^i(„+^_^) r(2/i + l) 



2i sin 7r(n + /i — (x;)r(l + n + /i — i^)r(l — n + /i + ti^) 
^ ^n+;.-^^^ _ ^^2m+2^i+i^^^^2^ + 1, 1 + n + /i + o;; 1 + n + /i - tl;; t). 

Proof. This follows from the generalisation of the Euler integral for the Gauss hyperge- 
ometric function and consideration of the consistent phases for the branch cuts linking 
the singular points, see pp. 91, section 17 " Verallgemeinerung der Eulersche Inte- 
grate". □ 

Remark 4.1. The first of these forms was given in |20|. 

For the description of H4.1^ in terms of the semi-classical form H3.2|l we have 

(4.4) m = 3, = {0,-1,-lA}, {p,}%i = {-^i-uJ,2uJl,2^i}, 

2V{Q) = -(m + Lj)t'\ W'{Q) = t-\ V{^t-^) = ^i^--r^ 



(4.5) 

and the coefficient functions are 

UN 



(4.6) 



Qn{z) = 



{N + 1 + IJL + UJ)Z 



rN 
rN+1 



{N + ll + Lj)t ^ 



r-jv 

rN+i 



{N + ii + Lj)z + {N+l + ii + uj)t ^ 



(4.7) 
(4.8) 

njv(2) = [1 + i(/i + a;)]z^ 

+ f (iV + 2 + /i + a;)(l - rM+ifN+i)^^ - + [1 + i(M + c^)]^ 

rjv+i Kjv+i t 

-[N+l{^^ + Lo)]t-' 
(4.9) 

^^^(^) =-Mm + '^)2' 



- CJi - - > 2 



- (iV + /i + tj)(l - rjv+ir-iv+i) 



^(/i-f Lj) 



1-Ht 



+ [N + 1 + i{ii + Lu)]r 



Our objective is to show that the average 111.311 can be evaluated via recurrence relations 
for the reflection coefficients. We begin with some preliminary lemmas. 

Lemma 4.2. The reflection coefficients for the weight {4 -If satisfy the homogeneous 
second-order difference equation 

(4.10) {N+l + ti + uj)trN+ifN ~{N + Lj)trMfN-i 

= [N + 1 + ^ + uj)fM+irN — (N -1+ + uj)fMrN-i- 
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Proof. The result of the above lemma, H4.1()^ can be found immediately from the general 
theory of Section 2, in many ways. By equating coefficients of z in the functional-difference 
equation 12.4711 using 14. 61 14. 714. 814. 911 . all are trivially satisfied except for the z coefficient, 
which is precisely 14.1UII . Similarly starting with (12.4811 and employing 14.6114.714.91 . one 
finds H4.f Alternatively one could start with either 12.49|l or l|2.5f |l and arrive at the 
same result □ 

Corollary 4.1. The sub-leading coefficients In, In satisfy the linear inhomogeneous equa- 
tion 

(4.1f ) (iV + ^ + ui)tlN - {N + ^ + uj)In = A'' [^(i - 1) + uj - ujt] kn- 

Proof. By substituting the general expression for the first difference of In, In using 12. IH 
in H4.1()|l one finds that it can be summed exactly to yield 

(4.12) (Ar + i + /i + tD)ti^ - (jV + i + /i + <.^)ii^ 

KN+1 I^N+1 

- {N + fi + u))t— + {N -\- fi -\- uj)— = ii{t ~ 1) + uj - Lut. 
Kn UN 

This can be summed once more to yield the stated result. □ 

Remark 4.2. One could alternatively proceed via the Freud approach |23| (see also |21|'l 
and consider the integral 



(4.13) / + + 



)J, + LU ^ ^ 2flt 



w{z)(I)n{z)4>n{z). 



1 + z 1 + tz 

Here we recognise the logarithmic derivative of the weight function in the integrand 

(4 14) 2^ 2^ 

^ ' w z ^l + ^^l + i^' 

and by evaluating the integral in the two ways we find a linear equation for In, namely 

Lemma 4.3. The sub-leading coefficients are related to the reflection coefficients by 
(4.15) 

In/hn + Un/i^n - N{t + 1) = ^ ~ "^^^^ [(iV + 1 + ^ + uj)trN+i + (A'' - 1 + m + •^)rN^\\ 

rN 

(4.16) = ^ ~ ''^''^ [(AT + 1 + ^ + a;)fjv+i + (A - 1 + M + (^)tfN-x\ 

rN 

Proof. The first relation follows from a comparison of the coefficients of z for Sliv(2:) given 
the two distinct expansions, the first by H,'i.8|l which reduces to l|4.8|l and the second by 
the specialisation of 1)3. 9|l . The second relation follows from identical arguments applied 
to VL*fq{z) or by employing l|4.10^ in the first relation. □ 

Remark 4.3. The first result appears in the Magnus derivation |35| for the generalised 
Jacobi weight, with 6\ = tt — <j>,92 = it, a = l-i,P = Wi,7 = —u!2. Then Eq. (14) of that 
work is precisely 14.151 . 
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Remark 4.4. The Magnus relation H4.15|l can also be found by employing the Freud 
method. In this one uses integration by parts on the integral 

(4.17) ^ 2^^"'(^ + + tz)'w'{z)4,N+i{z)j;A^, 

and in the term involving one employs (12.3511 for the derivative and 14.61 . 14.81 

for the coefficient functions. Equating this expression to a direct evaluation of the integral 
then yields (14.151 . 

Lemma 4.4. The sub-leading coefficient In can be expressed in terms of the reflection 
coefficients m the following ways 

(4.18) 24-^ = (iV + 1 + ^ + Lo)tC-^ ~ rN+ifN) + {N-l + ^l + lo)"^ 

KN TN rN 

-{N - 1 + ^ + iu)rNfN-i + iN + ti-uj)t + N-fi + iu 

(4.19) = (Af + 1 + At + ujf-^ + {N-l + fi + Cj)t{^^^ - rNVN^i) 

tn rN 

-{N + + Lu)trN+ifN + {N + fi- uj)t + N - fi + u! 
as well as analogous expressions for In- 

Proof. The first expression follows from a comparison of the z° coefficients for Qn{z) 
evaluated using both l|3.4^ and l|3.5^ . The second relation follows from an applying the 
same reasoning to Q*n(.z)- 

We will refer to the order of a system of coupled difference equations with two variables 
r„, fn say as q/p where g £ Z>o refers to the order of rn and p £ Z>o refers to the order 
of f„. 

Corollary 4.2. The reflection coefficients of the OPS for the weight \4. i| ) satisfy the 2/2 
order recurrence relations 

trNfN-i + riv-irjv ~t — l = L^fLlL [(jv + 1 + A^ + uj)trN+i + [N -1 + + u))rN-\\ 

rN 



(4.20) - [{N + At + uj)fN + (Af - 2 + At + Cj)tfN- 

rjv-i 

1 - rjvfjv 



[(iV + 1 + Ai + uj)fN+i + (iV - 1 + Ai + Cj)tfN-i] 
rN 

(4.21) - l-^^-i'^^-i [(TV + ^ + Co)trN + {N-2 + ^l + a;)rjv-2] 

rN-i 

and those specifying the solution for 11. 'X have the initial values 

(4.22) ro = fo = 1, ri = -w^^/wo, f-^ = -w^/wo, 
where the Toeplitz matrix elements are given in 

Proof. Solving (14.151 for the combination of In, In and differencing this, one arrives at 
(I4.2U1 . This however is of order 3/1 but by employing (I4.1U1 we can reduce the order in 
rjv of the recurrence to second order. The other member of the pair (14.211 is found in an 
identical manner starting with (14.161 . □ 
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Remark 4.5. The second-order difference H4.2()^ also follows immediately from equating 
the polynomials in z arising in the functional-difference l|2.42^ . after employing I|4.6I4.8^ . 
The other member of the pair, I4.21II . follows from the functional-difference 12.441 1. after 
using M. 714.911 . 

Remark 4.6. In their most general example Adler and van Moerbeke also considered this 
weight. In terms of their variables we should set P\ = P2 = 0,cii — t~^''^,d2 = t^^^, 
and without loss of generality 7" = 72 = 0. For the other parameters ^ = jj. — uj,^'i — 
2cJi, 72 = 2/^. There is a slight difference in the dependent variables due to the additional 
factor of t, so that we have the identification xm = { — l)'^t'^''^rM, J/jv ~ { — l)^t~'^^^fM 
and t;jv = 1 — rjvfjv. Generalising their working one finds that their Eq. (0.0.14) implies 

(4.23) -{N + l + ^ + Lu)xN+iyN + {N + 1 + fi + uj)xMyN+i 

+ {N - 1 + fi + uj)xNyN-i - {N - 1 + fi + uj)xN-iyN = 0. 

Now by transforming to our rjv, fjv and employing (12. IH one finds this is precisely 14.101 . 
which we showed is solved by 14.1111 . Their inhomogeneous Eq. (0.0.15) now takes the 
form 

(4.24) - 'UJV [(iV + 1 + + Lj)xN+iyN^l + N + fl + Lj] 

+ VN-1 [{N -2 + fl + Lu)xNyN-2 +N-l+^ + Uj] 
+ XNyN-l{xNyN~l + t^^^ + t^^'^) 

= -vi [(2 + + Lo)x2 + 1 + + oj] + a;i(a;i + t^^'^ + t~^'^). 
Upon recasting this into our variables and manipulating, it then becomes 
(4.25) 

trNTN-i + rNVN-i —t-1 "^NrN _|_ 2 + ^ -I- t^)fr]v+i + {N ~ 1 + fi + ci;)rjv_i] 

rjv 



1 - rjv-irjv- 



-[{N + fl + Uj)fN + {N -2+fJ. + Lu)tfN-2] 



rjv-1 

_ 1 - (1 - nfi) [{2 + + ui)tr2 + 1 + + Lj] + ri{tri -t-1) 

rNfN-i 

However using the identity 

C2Fi{a,h\c;x) = [c+ {1 + b- a)x]2Fi{a,b + l; c+1; x) - ^^-^(1 + c- a)x2Fi(a,b + 2; c+2; x) 

c + 1 

we note that the right-hand side is identically zero for the initial conditions l)4.22|l and the 
recurrence is not genuinely inhomogeneous, thus yielding our first relation above, H4.2()^ . 

We seek recurrences for rjv, fjv which are of the form of the discrete Painleve system 
ifLSji . ltL9|l . For this purpose a number of distinct forms of the former will be presented. 
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Proposition 4.1. The reflection coefficients satisfy a system of a 2/0 order recurrence 
relation 



(4.26) 

|(1 - r-ivfiv) [{N + 1 + ^, + Lj){N + ^. + uj)trN+i - (N + ti + Lu){N - 1 + fi + uj)rN-i, 

+ N{N + 2iJi){t~l)rN] 
X I (1 - rNrN) [{N + 1 + ti + uj){N + + uj)trN+i - {N + + uj)(N - 1 + fj. + uj)rN- 
+ (TV + 2^)(iV + 2^ + 2ui){t - l)rjv} 
= -{2N + 2At + 2wift{l - rNfN) 
X [(iV + 1 + /i + uo)rN+i + (iV + A* + ^^Yn] [(iV + a* + ^^Vn + (iV - 1 + a* + t^)'-iv-i] 



and a 0/2 order recurrence relation which is just \4.2b\l with the replacements u ^ u and 



Proof. Consider first tlie specialisation of l|8.44|l to our weigiit at liand at the singular 
point z — —1, and we have 



(4.27) 1-^ - Nt-' -(iV + l + A. + a;)^^ + - t~')] 
( UN Kir rN J 



+ 



2 

KjV-l 



Ar + /i + a; + 



{N-l + fi + iu) rjv,i 



rN 



(^ + ^ + ") +(iV + l + M + ^)'^ 



2 / i-1 



rN 

2 



by using I|4.6I4.8^ . Similarly l|3.44^ evaluated at z = —1/t yields 



(4.28) |iil-7V-(iV + l + A. + ^)^^ + f,{t-' - 1)1 
N + ^ + Ld + {N-l + ^ + Lj 



2 

AT 



, »"iV-l 



N + ll + LU+{N+l+fI + Lj] '"'^ + ' 



t - 1 



rN 

\ 2 



The first relation follow by eliminating In between 14.2711 and 14.281 1. whereas the second 
follows from an identical analysis to that employed in the proof of Proposition 14.41 but 
starting with the bilinear identity 13.45t . □ 
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Proposition 4.2. The reflection coejficients also satisfy a system of 1/1 order recurrence 
relations the first of which is 



(4.29) 

I - (Af + A* + - rNfN)t [(iV + 1 + ^ + Lu)rN+ifN + (iV - 1 + ^ + uj)rNfN~i] 
+ 2{N + M + i-jfrjffM -{N + fj. + ujf{t + l)r-jvfjv - 2{N + ^ + uj)uj{t - l)rNfN 

+ {ti-cj){ti + cj){t-l)} 

X I - (Af + + - rjvfiv)t [(TV + 1 + ^ + u!)rN+ifN + (N - 1 + + Ld)rNfN-i] 
+ 2{N + A* + i^fr%f% -{N + fi + ujf{t + l)rjvfjv + 2(jV + ^ + uj)uj{t - l)rNfN 
+ {p - uj){fi + Lu){t - 1)^ = - [2{N + + oj)rMfN + (jJ-ujf{l- rNfN) 
X [(iV + 1 + /i + uj)trN+i + (TV + ^ + uj)rN] [(A^ + A* + + (AT - 1 + a* + uj)tfN^i\ 



and the second is obtained from {4 ■ 29^ 1 with the replacements ui ^ ijj and t^^^^Vj 
t^^l\.. 



Proof. The specialisation of 13.461 to tlie weiglit 14.11 evaluated at the singular point 

2 = — 1 is 



(4.30) 1-^ - Nt-' + iN + ^, + Lj)t-'^ + a;i(l - t-')Y 

2 

+ [{N + 1 + 11 + iJj)rN+i + {N + ^1 + a;)t" Vat] 

[(Af - 1 + M + Oj)fN-i + {N + ^i + a;)t~Vjv] = ujI (^-^) 



by using Il4.(il4.8t . Similarly Il3.4()t evaluated at z = — 1//; yields 



(4.31) {hL^N + {N + ^, + u)^ + ^,{t-^-l)y 

I KN H at > 



^2 ^ 

+ [{N +1 + ^. + Co)rN+i + {N + 11 + Lu)rN] 



X 



[(Af - 1 + ^ + uj)fN-i + {N + ^ + cj)fjv] = 



t-1 



Again eliminating In between these two equations yields the recurrence relation 14.291 . 
The second follows in the same way starting with 13.471 . □ 
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Proposition 4.3. The reflection coefficients satisfy an alternative system of 1/1 order 
recurrence relations the first of which is 

(4.32) [(iV + l + ^ + tD)(Af + /i + iI;)triv+ifiv 

-(iV+l+^ + Cj)(iV + /^ + Lj)fN+irN + + n){t - 1)] 

X [(iV + 1 + ^ + tj)(iV + + uj)trN+irN 
-(iV+l + At + '^)(^f + M + ^)fN+irN + (c^ - m)(^^ + At)(* - 1)] 

— (id ~ LO)^ 

X [(AT + 1 + + Oj)trN+i + (Af + A* + t^)riv] [(iV + 1 + ft + (.^)riv+i + (W + ^ + a;)ffiv] 
and £/ie second is again obtained from \4-.3'^ with the replacements ui uj and t^^^^rj ^ 

Proof. The specialisation of H3.48^ to the weight l|4.1^ evaluated at the singular point 

2 = — 1 is 

(4.33) ( -^^^ + (iV + M + u;)fN+irN +u}i + (p~ iuj2)tY 

L KJV + I J 

= [{N+l + fi + uj)trN+i + (iV + M + ^Vn] 

X [(iV + 1 + ^ + Uj)fN+l + (iV + A* + '^)tfN] +LJi{t- if, 

by using 14.614.81 . Similarly 13.481 evaluated ai z = —1/t yields 

(4.34) (-^^^ + (Af + ^ + a;)fjv+irjv + t<:' + /i4^ 

= [{N+l + fi + iu)rM+i + (AT + ^ + cj)rjv] 

X [(Af + 1 + ^ + Lu)fN+l + {N + fi + tj)fjv] + fi'^it- if, 

Again eliminating Zjv+i between these two equations yields the recurrence relation 14.321 . 
The second follows in the same way starting with 13.491 . □ 

Remark 4.7. Note that the recurrence system 14.261 and its partner is quadratic in r]v+i, 
TAT-i and fM+i, rN-i, the system 14.291 and its partner is also quadratic in rM+i,rN-i 
and f]v+i, rjv_i, and likewise 14.321 is quadratic in rjv+i, fjv+i. This renders them less 
useful in practical iterations than the higher order systems that are linear in the highest 
difference. By raising the order of one of the variables by one we can obtain a recurrence 
linear in the highest difference. 

Corollary 4.3. The reflection coefficients satisfy a system of a 2/1 order recurrence 
relation 

(4.35) {N + 1 + H + - uj)t{l - rjvfjv)r]v+i 

+ {N -1 + H + uj)[2{N + fi + Lu)rNfN +1^0- Lj]rN-i 
- {N -1 + n + Ld){2N + 2/i + 2L0i)tr%fN-i 

+ [{uj - L0)N{t+ 1) - (2^4 + 2uJl)[ll{l -t) +LOt- Ld]]rn = 
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and a 1/2 order recurrence relation which is again obtained from \4..S5\l with the replace- 
ments U! <-i oj and t^^^^rj ^ t^^^^fj. 

Proof. The solutions for the sub-leading coefficient In, In that arise from the simultaneous 
solution of 14. 3014. 3111 and 14.3314.341 respectively are given by 

(4.36) 

=\{N + ^i + uj)t{l - r-ivfiv) [(iV + 1 + ii + iu)rN+if=N + {N - 1 + + iu)rNf=N-i] 

KN ^ 

+ {N + ^. + ui)[N{t+l) - ^(1 -t)- ujt^iL]rNfN + (tJ + /i)[/i(l -t)+uot- 
[2(A'' + + uj)rNfN + u) — u] 

(4.37) 

^{{N + ^Ji + u) [(TV - 1+ ^ji + uj)trNfN-i - {N -1 + ^i + uj)fNrN~i\ 

UN L 



+ {lo -\- — t) -\- ujt — uj]^ I (lo — uj), if tj 7^ tj, 



and the corresponding expression for In/i^n under the above replacements. Equating 
these two forms then leads to 14.351 . □ 



The systems of recurrences that we have found are in fact equivalent to the discrete 

4 ^ ^5 



Painleve equation associated with the degeneration of the rational surface D^^ — > Dp' 



and we give our first demonstration of this fact here. 

Proposition 4.4. The N -recurrence for the reflection coefficients of the orthogonal poly- 
nomial system with the weight \4- if is governed by either of two systems of coupled first 
order discrete Painleve equations 1 j.^) . U.ytl . This first is 

,ifN + N)ifN-^N + 2^i) 



(4.38) QN+igN = t- 



/jv(/Ar - 2l^i) 



(4.39) /jv + fN-i — 2.L01 H 1 , 

gjv — 1 9N — t 

subject to the initial conditions 

(4.40) s^-^ 'l^lnl'lt' ' = 

^i + ct;+(l + /x + Ljjtri 
The transformations relating these variables to the reflection coefficients are given by 

N -1 + fi + uj + {N + fj. + uj)-^'^ 



(4.41) gN = t 

N-l + fi + Lu + (N + n + uj)t 

riv-i 

(4.42) fN 



t^^N-{N+l + ^i + (I;)(l - rNfN)t''''^^ 



KN rN 



1 - 1 

The second system is 

ifN+N){fN+N + 2uji) 



(4.43) gN+iQN = t ^- 



/iv(/jv - 2fj,) 



(4.44) fN + fN-i =2fi+ _ ^ + ^ 1^—-^ , 

gN - l qn - t ^ 
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subject to the initial conditions 

H + LJ + {l+fl + ljj)t^^fl 



(4.45) 



gi = 



/o = 0. 



fi + uj + {l + fi + uj)ri 
The transformations relating these variables to the reflection coefficients are given by 



N-l + ^ + Lu + {N + fi + uj)t' 



-1 TN 



(4.46) 



(4.47) 



QN 



In 



N -1+ fl + LU + {N + fl + Lu)-^ 

rjv-i 



-thL + Nt + {N -l + fi + Lu){l~ rMfnV-^ 

KN TN 



Proof. Consolidating each of 14.2711 and 14.281 into two terms and taking their ratio then 
leads to (14.3811 after utilising the definitions 14.4114.4211 . The second member of the recur- 
rence system 14.391 follows from the relation 

(4.48) ln+L + hL 

K.N+1 UN 

{N + 2 + fi + uj)(l- rjv+ifjv+i) + (iV + M + io)t'^^^ -(N + l + fi + iu)rN+if=N 



vn+i 



rN+i 



2lui - 2^lt-^ + (AT + 1 + ^ + (I;)(l + t~^), 



which results from a combination of 14.181 and 12.111 . and the definition 14.421 . All 
the results for the second system follow by applying identical reasoning starting with 

lEZSi. □ 

Remark 4.8. Generalised hypergeometric function evaluations were given in j22| in the 
special case ^ = 0. In terms of our unitary group average one such evaluation reads 



(4.49) mzr-^ii+zif^K^+t^i) 



Urn 



2n\ 

/ U(N) 



3 = 



(i + M + ^ + i)r(i-M + ^ + i) 



X 2-P'i^'(~2M; ~M - ^ - M + ■ • • ,iiv)|ti=. 

subject to 5R(tJi) > — V2 and \t\ < 1. Similarly, for the reflection coefficients we have 
(4.50) 

-/i — tj;7V+l — /i-|-a;;ti,...,tjv) 



^/ .-.N (M + ^)iv 2J'i^'(-2m. 1 



F^'> (-2^, -fi-u;N - n + uj;ti,...,tN) 



(4.51) 



{—fJ, + Lj)is[ j-P"!^"* —1 — fi — ui; N — 1 — fi + uj;ti 



, *iv) 



{l+fl + Lj)l 



2Fi^\-2fi,-fi-uj;N - fi + ij-ti 



1 1 — . . .— t w— i 



The analog of the Euler identity for this function is 
(4.52) 

2Fi^\~2n,-fi-uj-N - ^i + Lu;ti, . . . ,t]v)|ti=...=t„=i 



n 



Tjj + 2ii + 2cji)r{j - fi + cj) 

r{j + 2oji)r{j + fi + Lu) ' 
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when 5R(/i + lji) > — V^i + lu) > —1, thus implying Gamma function evaluations in 

the special case t = 1. 

Remark 4.9. The degeneration of the Pvi system to the Pv system is facilitated by the 
replacements aj + /^i-^!/, a) — /ih^/i, in^ and then taking the limit /i oo. The 

coefficients of the orthogonal polynomials rN,lN remain of 0(1) in this limit. Then we 
see the explicit degeneration of the following equations - 1)4.11^ — > Eq. (4.23) '21', the 
recurrence relations l)4.29|l —^ Eq. (4.60)|l21j, l|4..'-{5^ Eq. (4.9)121] and its conjugate to 
Eq. (4.10)[22 modulo the identity Eq. (4.5)123, a-^d the hypergeometric functions H4.49^ 
^ Eq. (4.24)121], Eq. (4.26)|2I|, and llI3Tt ^ Eq. (4.27)|2I]. 

Remark 4.10. Two simple cases exist for the special values of the argument t — 0,1. In 
the first case, t = 0, we have 

(4.53) r^^i-lf-^t±^^, = (-^ + -)^ 



(1— /1 + Lj)jv' {l + IJ. + Uj)l 

whereas for t = 1 (and ^ is irrelevant) we have 

(4.55) rj, = {-lf-^^^±^^, fr, = {-lf + '^^^ 



(4.56) In = 



(l + /i + a))jv' (l + n + uj) 



{N + ^I + LO)' 



There is a specialisation of the generalised Jacobi weights leading to a formulation of 
the orthogonal polynomial system in terms of real variables, and a simple phase factor 
appearing in the reflection coefficients. This occurs when £ R, tJ2 = and \t\ = 1 

and is a special case of hermitian Toeplitz matrix elements. In such a situation r„ is no 
longer independent of r„ (it is the complex conjugate of r,i) and the coupled systems of 
difference equations reduce to a single equation. 

Corollary 4.4. When L02 ~ 0, t £ T , fi,uji £ M. with jj, + uii ^ Z<o and fi = tri then the 
reflection coefficients are products of a real coefficient a;„ £ R and a phase factor so that 

r — f-^l'^T r — t^l'^T 

Proof. Setting lj2 = in H4.1()|l we note this can be rearranged as 



+ (n - 1 + ^ - 



= 0. 



Given that r-„ — t"r„, r,i_i = t" ^r„_i we use the above equality to show 
(4.58) f„+i=r+V„+i, 

and by induction on n the statement fn = t"r„ must be true n > as it holds for n — 0,1. 
The corollary then follows. □ 
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5. The t-function Theory for Pvi 

In the Okamoto theory for Pvi the Hamiltonian function which governs the evolution 
of {q, p; H, t} through the system 11.711 is 

(5.1) K:=t{t~l)H 

= q{q - 1)(9 - t)p^ - [a4(g - l)(g -t) + azq{q - t) + (qq - l)q{q - l)]p 

+ Q2(ai + Ci2){q - t), 

with parameters ao, ai, Q2, Qs, Q4 G C subject to the constraint qq + ai + 2ot2 + Q3 + 
0:4 = 1. In the introduction it was remarked that the unitary group average 11.31 was 
shown to be a r-function for the sixth Painleve system |22| and this can be achieved via 
two distinct methods. In the first method I2UI the connection with the Pvi r-function 
was established for an average with respect to the Cauchy unitary ensemble (see Eqs. 
(1.12,1.19,3.20,3.28,3.30) of |22| ') and using the stereographic projection this average was 
related to 11.311 with the parameters 

(5.2) (qo, Qi, Q2, Q3, Q4) = + 1 + 2a;i, iV + 2^i, -iV, -fi-uj,-fi-Lj 

The appropriate sequence of the Hamiltonian variables {qn,Pn, Hn,T,i},i=o,i,... in which 
A'' is only incremented is generated by a shift operator L^^ — riSoS\S2S3S4S2 in terms of 
the reflection operators and Dynkin diagram automorphisms of the extended afRne Weyl 
group WaiO'l^^) (see [221). It has the action L^^ : ao 1— > ao + 1, ai > ai + 1, 02 1— > Q;2 — 1. 
For such a sequence we have the following result. 

Lemma 5.1 f|2()|.|22|'l. The sequence of auxiliary variables {gn, fn}n=o,i,... defined by 

(5.3) gn ■■= , 

(5.4) /„ — qniqn ~ l)pn + {1 - 02 ~ a4,)iqn - 1) - asgn - Qo ^"^'^ — 

qn — t 

generated by the shift operator Lq^ satisfies the discrete Painleve equations associated with 
the degeneration of the rational surface D'^^'^ — » D'^'^ 

frc-\ _ * (/n + 1 - Q!2)(/n + 1 - a2 - 04) 

(O-O) gn + ign 



(5.6) f„ + /„_i = -Q3 + 



t~l fnifn+as) 

ai aot 



gn-l t{gn -1) - gn 

Applying this result to the average ifTTT^ impl ies a recurrence scheme to compute the 
latter. 

Proposition 5.1 ([201). Let {(7jv}jv=o,i,..., {/jv}jv=o,i,... satisfy the discrete Painleve 
pled difference equations associated with the degeneration of the rational surfa 

t (/jv+iV + l)(/iv + iV + l + M + ^) 



cou- 

) _^ 

4 ~^ 



(5.7) gN+igN = 



(5.8) /jv + /iv-i =M + '^ + 



t-1 fNifN-fi-Uj) 

N + 2(1 {N +1 + 2uji)t 



gN - 1 t{gN - 1) - 5JV 
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where t = 1/(1 — e"*) subject to the initial conditions 

go = fo = {l+f^ + 0)iqo - 1) + {li + Lu)qo - (2a;i + l)S!liS!LzJl 

qo - 1 qo -t 

with 

(5.9) ^" = Kl + ^i^°S^"''^^(^'') 

Define {qN ,Pn}n=o,i,... by 
(5.10) 

n - 

In — -, 

gN -1 

(5.11) 

_ (ffiv - 1)^ 
PiV — /jv 

- (iV + 1 + M + Lo)^^^-^ - (m + c^)(5iv - 1) + (iV + 1 + 2c^i)- ' ^ 



gjv t+(l-t)g]v' 

Then with ro(e''*) = 1 and Tx{e''*') = «;o(e''*) as given by 14.^14.3^ , {Tn}n^2.3,... «s 
specified by the recurrence 

(5.12) - (TV + M + ^)(iV + M + ^) 

= g]v((7iv - l)pw + (2^i + 2uJi)qNPN - (m + i^)PiV - A''(A'' + 2^ + 2iUi). 

In the second method the connection with the Pvi r-function was estabhshed for an 
average with respect to the Jacobi unitary ensemble (see Eqs. (1.21,3.7,3.28,3.31,3.32) of 
|22|'l and using the projection ( — 1,1) T under the condition ^ = this average was 
related to 11.311 with the parameters 

(5.13) (ao, Qi, a2, as, 04) = (^l — fi — uj, N + 2fi, —N, -fi — lj, N + 2uji^. 

Sequences of the Hamiltonian variables {qn,Pn, H„,t„}„=o,i.... are now generated by the 
shift operator L]~/ = r3SiS4S2SoS3S2. It has the action L^^ : qi qi + 1, 02 1-^ 02 — 
1, 04 ^ Q4 + 1. Using the methods of |22) we have the following resuH. 

Lemma 5.2. The sequence of auxiliary variables {(?n, /n}n=o,i,... defined by 

(5.14) g„:^-^"~' 



(5.15) /„ 



. - 1 
1 



1 - t 



{qn - t){qn - l)p„ 

(<?« - i){qu - 1) 



+ (1 - ao - a2)(gn - 1) - a3(g„ - t) - a4- 



qn 

generated by the shift operator L~[l satisfies the discrete Painleve equations associated with 
the degeneration of the rational surface D'^^'^ D'^'^ 

/c-i«\ , (/n + 1 - a2)(/„ + 1 - ao - 02) 

(5.16) Sn + lffn = t ^ 

/n(/n + 03) 

(5.17) /„ + = -as + + 
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Proof. Using the action of the fundamental reflections and Dynkin diagram automor- 
phisms given in Table 1 of |22| we compute the action of L^/ on q and write it in the 
following way, 

jS. — M| ^ = t[q{q- l){q - t)p+{ai +02)9^ - ((qo + ai + 02)* - ao - a4)q - a4t] 

X [q{q - l)(g - t)p + (ai + 02)9^ - ((ai + a2)t - 04)17 - ctii] 

[q(q — l)(g — t)p + (ai + 02)?^ — {—ait + ai + a2)q ~ ait] 

[q{q - l)(q - t)p + (fti + a2)q^ - {-{as + 04)* + ai + Q2 + Qa)? - cat], 

where q := q„,q := qn+i- From the definitions 15.1415.151 this result can be readily recast 
as IS.lfcil . The second 15.1711 follows from a computation for /„ + /n-i using the shift 
operator L14. □ 

Remark 5.1. The two systems of recurrences H5.5I5.6^ and H5. 1615.17(1 are related by an 
element of the 54 subgroup of the Wa{F4) transformations, namely the generator |48) . 
This has the action 

(5.18) : ao ^ a4,t J~[''i ^ ^ ~ ^^"P' 

and when applying these transformations to 15. 31 . 15.41 . 15.51 . 15.6(1 we recover 15.141 . 
lIKTHll . lEHnj, iICTtIi respectively. 



Proposition 5.2. Let {givjiv^o.i,..., {/]v}]v=o,i,... satisfy the discrete Painleve coupled 

the degeneration of the rational 
{fN+N +l){fN+N + ti + uj) 



difference equations associated with the degeneration of the rational surface D^^'' —» Dg^' 



(5.19) gN+igN = t- 



(5.20) /at + /iv-i = M + + 



./iv(/iv - M - ^) 

TV + 2^ {N + 2iUi)t 



ffjv - 1 Qn -t 
where t — e"^ subject to the initial conditions 



90 - i „ 1 
50 = 7, fo - 



go-1 1-t 
with 

_.d_ 

Ldl d(f> 



{^l + UJ){qo-l) + {^l + ^){qo~t)-2uJ^^'' 



J —log e^'"l'Ti{e^'l') 



(5.21) qo = 



d 



^ /i + (.j + i-5^1oge*'"^ri(e"*) 
d<p 

Define {gjv,pjv}jv=o,i,... terms of {f^ , Qn} n=o.i,... by 

gN -t 



(5.22) = 



9N -1 



(5.23) piv = jrhv- I'^^N - l)fN - (m + ^)9N + {N + 2cu^) ^^ -N-^.-u^. 

(1 — t)gN L 9N — t J 

r/ien m£/i ro(e''*) = 1 and Ti{e'''') = wo(e''^) as given by (T^T^, {r]v}]v=2,3,... is 
specified by the recurrence 

(5.24) -{N + ^ + ^){N + ^ + u;)^^^^±g^ 

= giv(gjv - lfp% + [{2fi - N)qN + N + 2oJi]{qN - 1)pn - 2fiNqN - iV(A'' + 2uji). 
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Proof. Let Y„ ■- L^^Kn - K„ ^ A"„+i - /C^ and from Table 1 of |^ we have 
{t-l)q„{, -.N , , 1 , (1 - ao - Q2)(ai + ^2 + as) 



(g„ - l)p,i + (oi + a2)gn + "4 
Now consider 

and this latter difference, upon again consulting Table 1 of 22 , turns out to be 

i(t — 1)^ log [qn{qn — + [(ai + 2a2)gn + «4](gn - l)pn + Q2[(ai + a2)qn + 04]] . 

After integrating both expressions and introducing an integration constant Eq. 15.241 
follows. □ 

We now seek to relate the results of the r-function approach to the theory developed 
for the orthogonal polynomials on the unit circle with semi-classical weights as given in 
the previous section. However we will only discuss the scheme given in Proposition 15 . 1 1 as 
this is the simplest. 

Proposition 5.3. The transformations linking the Hamiltonian variables qN,PN m Propo- 
sition \5J\ to the reflection coefficients rN,fN for the system of orthogonal polynomials with 
the weight {4- -?| ) are given implicitly by 

{N + ^1 + Ld)rNfN 1 

gjvpjv + ^ + (.tJ = 



(TV + ^ + uj)rNrN — /J, + UJ qN — i 



(5.25) 
(5.26) 

(5.27) 



{N + 2a;i)(gjv - 1) - t— + Nt + {N + 1 + fi + - rjvriv)t''''^+' 



= {N + fj, + u;)[{N + + uj)rNfN ~ fi + i^u] 

X -, '-^ 

{N + 2iJi)qN +t—-Nt-{N-l+n + i:u)(l- rjvfiv)*^^ 
KJV rN 

[qN — 1)pjv + h + uj = [N + n + io)[{N + fi + Lu)rNrN — fi + uj] 

qN 

X -. 

(TV + 2iUi){qN -l)~t—+Nt+{N + l + fi + ij){l- rjvfjv)*^^ 

KN rN 

_ {N + fi + uj)rNfN 1 



{N + fj, + Lo)rNrN — fi + Lu qN 



(5.28) 



{N + 2LJi)gjv +t— - Nt-{N -1 + fi + Cj){l^ rjvfjv) 



Proof. We require in addition to the primary shift operator L^^ generating the A'^ 1-^ A'^ + l 
sequence another operator which has the action iuj2 1-^ iL^2 — 1- This is the secondary shift 
operator T^^ — r\SiS2SoSiS2Si and has the action T^^ : aa — > Q3 + 1, 04 — > Q4 — 1. From 
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Table 1 of |22) we compute the actions of Tg^^ , T34 on the Hamiltonian to be 

r34^ ■ K„ - Kn = ~ q„{qn - l)Pn 

+ (qo + a4 - l)(gn - 1) - (Q2 + a3)(ai + 02 + Q3)- 



(g„ - l)p„ - 03 



134 ■ K„ - K„ = - q„{q„ — l)p„ 



+ (qo + 03 — l)gn - (Q2 + Q4)(ai + 02 + 04) 



^nPn — Q4 



However 



. - = t(f - 1)^ log § = t(t - 1)^ logrn 



and we employ the results of Corollary 13.21 and the evaluation of the coefficient functions 
in 114. 814. 91 to arrive at 

TN UN rjv 

(t - 1)— = — ^ + iV+(iV-l + ^ + io){l - rjvf^)^^. 
rjv UN Tn 

In addition we note that after recalling 1)2. 9|l . 1)5.12^ factorises into 

{N + ^ + u)){N + + ui)rMfN = [qjvpjv + ^ + ti)] [(giv - l)p]v + ^ + c^;] . 
The stated results, 115. 2515. 2811 . then follow. □ 

6. Applications to Physical Models 

6.1. Random Matrix Averages. A specialisation of the above results with great inter- 
est in the application of random matrices |32| is the quantity 

N 

(6.1) FS''''{u;f^):^(ll\u + z,f^) 

\ f-J-^ I CUEjv 

This has the interpretation as the average of the 2/i-th power of the absolute value of the 
characteristic polynomial for the CUE. In the case |u| = 1 16.11 is independent of u and 
has the well-known (see e.g. |S]) Gamma function evaluation 

AT AT JV-1 



^ ^ ' /cUE„L=e..* \\}^^ " /CUE„ r2(i+l+At) 

when K(^) > — V^- For |it| < 1 we see by an appropriate change of variables that 

JV JV 

(6.3) (ni"+^'i'")._ =(n(i+i"^o"(i+iM)'' 



CUEjv \ I CUEjv 

(6.4) = 2-Fi'^'(-M, ~fJ-;N-ti, . . . ,tN)\ti=...^t,^^\u\^, 

where the second equality follows from 14.4911 . For |u| > 1 we can use the simple functional 
equation 

JV JV 

|2(j\ _| |2;iJV/ TT I 1 |2n\ 



(6.5) (ni-+-'i''^) =i«r^(ni-+ . . 

I CUEjv CVE.1, 

to relate this case back to the case \u\ < 1 
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The weight in the first equality of l|6.4^ is a special case of H4.1|l . In terms of the 
parameters of the form H4.1^ we observe that = 0, 2fi ^ n, uj — ui = /i/2, i.e. L02 = 
and t = I It p. The trigonometric moments are 

(6.6) w-n = iw^^,^/^ — -2Fi{-ii,-ii + n;n + l;\uf) n G Z>o 

nil (n + 1 — n) 

(6.7) Wn = \uf'^W-„ n G Z>o. 

The results of Section 0] then allow 1)6. l|l to be computed by a recurrence involving the 
corresponding reflection coefficients. 

Corollary 6.1. The general moments of the characteristic polynomial | det(ti + U)\ for 
arbitrary exponent 2fi with respect to the finite CUE ensemble U £ U{N) of rank N is 
given by the system of recurrences 

pCUE pCUE 

,„ r,^ ^N+l ^N-l , I |2JV 2 

(6-8) (j.cuE)2 =1-1^1 ^JV' 

with initial values 

(6.9) ^0^^™ = !, Ff'™ = 2Fi(-M,-M;l;hl'), 

and the recurrence relation for the reflection coefficient rjv 
(6.10) 



1 - Inl'^'^r^ 

2\uf''rNrN^i - luf - 1 = ^ [(iV + 1 + fi)\ufrN+i + {N - 1 + f^)rN-i] 

rjv 

_i_ju\ [(;v + ^)|^|2^^^(^_2 + M)r^_2] , 

rjv-1 

subject to the initial values 

, 2i^i(-M,-M + l;2;|^|^) 
6.11) ro = 1, ri = — I ■ 

2-Fi(-M, ) 

Proof. From either 114.11)11 . 14.1111 or 14.351 1 and the fact that fi — \u\^r\ we can repeat 
the arguments of Corollarv 14.41 to deduce that fjv — jup^rjv for A'^ > 0. The recurrence 
relation follows simply from the specialisation of 14.201 and the initial conditions from the 
= 1 case. □ 

Another spectral statistic of fundamental importance in random matrix theory is the 
gap probability for the circular unitary ensembles, and this is the specialisation whereby 
fi = uj — ijj = 0, |f| = 1 so the angle (j> G [0, 27r), whilst ^ G C is general (one is mainly 
interested in an open neighbourhood of ^ = 1). The generating function for the probability 
of finding exactly k eigenvalues z = e'" within the sector of the unit circle 9 £ {n — (f),TT] 
is denoted by E'^^{{0, 4>);^) and has the definition 
(6.12) 

1 / fiv fir \ / PTV fir \ 

i5^^^((0»;O-7^( / -W -d n le'^-e'"]", 

where the normalisation Cjv = (27r)^A''!. It is well known that Toeplitz elements with 
such a symbol have the form 

(6.13) =5n,0 + 7r^(-l)" 
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which is easily recovered from the general expression H4.2|l . A recurrence scheme for the 
generating function H6.12|l involving particular examples of the coupled discrete Painleve 
equations 11.81 . 11.91 has been presented in |2U| . Here we use recurrences found herein for 
^■jv, tn, together with the fact that for 16.121 one has rjv — t~^fjv, to replace the role of 
the coupled recurrences from |2()| by a single recurrence. 

Corollary 6.2. The generating function for the probability of finding exactly k eigenvalues 
z = e'® from the ensemble of random N x N unitary matrices within the sector of the unit 
circle 9 G (tt — tt] is given by the following system of recurrences in the rank of the 
ensemble N, 

TpCVE 7T.CUE 

(6-14) (^CUE)2 =^-^N, 

where the initial values are 

(6.15) -Bo = 1, El = 1 - — 

and the auxiliary variables xn are determined by the quasi-linear third order recurrence 
relation 

(6.16) 2XNXN-1 - 2 cos ^ = }_^zlk [(iv + 1)xn+i + {N ~ 1)xn-i] 

2 XN 



_ 1 - X%_i 
XN-1 

or the quadratic second order recurrence relation 



[NxN + iN~- 2)xN- 



2 



(6.17) (1 - xlf [{N + ifxl+r + (TV - lfxl_^] + 2{N^ ~ 1)(1 - x%)xn+iXn-i 

2 (t> 2 1 n 
2^ = 0, 

along with the initial values 

. 



+ 4iVcos |a;jv(l - xif) [(iV + l)a;jv+i + (iV - 1)xn-i] + ^N^xj, 



(6.18) x-i =0, xo = 1, xi 



2n^ 



Proof. The first recurrence relation follows directly from the general recurrence ijOnjl and 
CoroUarv 14.41 whilst the second follows from 14.2911 . □ 

6.2. 2-D Ising Model. It has been known for some time that the diagonal spin cor- 
relations in the square lattice Ising model could be evaluated in terms of the Painleve 
sixth transcendent |29) and in this work a coupled system of difference equations involv- 
ing eleven variables were given. Here we give what we consider to be the simplest set of 
recurrence relations for these correlations as a special case of the general theory above. 
The diagonal correlation functions are given by |42| 



det(ai_j(fc))i<ij<]v if fc > 1 or T < Tc, 
det(ai_j(fc))i<ij<]v if fc < 1 or r > Tc, 



(6.19) (o-Q.oo-iv.iv) = 

where 

1 f I i-^-i — 1 1 /■ 

dec 



1 - fcc 
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and the argument is defined k — sinh^(2 J/fcsT), J being the coupling strength and T 
the temperature. The weight appearing in the Toephtz determinant form of the low 
temperature correlation is 

(6.21) w{z) = Cz'/^ll + z\'^^^{l + k-'^zf''^ = Cz^''^(\ + + fc-2^)'/^ 

which is a special case of 14. H with ^ = 0, /i = 1/4, oji — —1/4,0^2 — i/2 and t — Xjk? . 
Here a„ = ( — l)"fe~"w_„(l/fc^). In the high temperature regime the exponents are 
sign reversed and t = k^, so a„ — ( — l)"fc"w_„(fc^). So these cases form an interesting 
example where UJ2 7^ 0, so both rjv,fjv are distinct and independent in contrast to the 
random matrix and quantum many-body cases. As is also well known the Toeplitz matrix 
elements in the low temperature regime are given by 

(6.22) n^_„ = t}LIh + yi}E^,F^{-y2,n + y2-,n + l-,k-^), n > 0, 

TV 1 (n + 1) 

(6.23) = (-l)"+^fc-^" r(n - V^)r(y2) ^^ _ y^.^^ 

TT 1 (n + 1) 

whilst those in the high temperature regime are 

(6.24) ^ + t^;r^ i-.(V.. + V.n + 2; n > 0, 

TT 1 (n + 2) 

(6.25) ^'n^ '^^""yy^'^'^ 2Fi(-V2,^-V2;n;fc^), n > 0. 

TTK 1 (71) 

Corollary 6.3. The diagonal correlation function for the Ising model valid in both the 
low and high temperature phases (with fc 1/fc in the latter case) is determined by 

,r-nc!\ (cr0,0O"JV + l,]V+l){o"0.0O"JV-l,JV-l) , 

(6.26) -3 = l-rjvr-iv, 

(cro,oo"jv,jv; 

along with the quasi-linear 2/1 

(6.27) (2iV + 3)fc"'(l - rivfiv)rjv+i + 2N [k'^ + 1 - (2iV - l)fc- Vjvf jv-i] rjv 

+ (2iV - 3) [(27V - l)rjvfjv + 1] r-jv-i = 0, 

and 1/2 recurrence relations 

(6.28) (2iV + 1)(1 - rNfN)fN+i + 2N [{2N - 3)fjvrjv-i + k'^ + l] fjv 

+ {2N - l)fc"^ [-(27V + l)rjvfjv + 1] fjv-i = 0, 

subject to initial conditions for the low temperature regime 

, , 2-fc2 fc2-lK(fc-i) ^ fc2-lK(fc-^) 

(6.29) ro = l, fo = 1, ri = — — + — — ) ' , fi = -1 ' ^ '- 



3 3 E(A;-i)' fe2 E(fc-i)' 

or to t/ie initial conditions for the high temperature regime given by 
(6.30) 

1/2 E{k) \ _ k^E{k) 

ro = 1, ro = 1, ri^-{ — ~ — —— n = 



3ifc2 (A:2 - 1)K(A:) + E(fc) j ' (fc2 - l)K(fc) +E(fc)' 

where K(fe), E(fc) are t/ie complete elliptic integrals of the first and second kind respectively. 

Proof. 116.2716.28^ follow from H4.35^ and its "conjugate" upon the specialisation to the 
parameters above. The initial conditions follow from explicit evaluation of the Toeplitz 
determinants. □ 
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The correlation function and reflection coefficients have particularly simple, yet general 
forms when expressed in terms of generalised hypergeometric functions. 

Corollary 6.4. In the low temperature phase the diagonal correlation function is given 
by 

(6.3f) (o-o.oo-iv.jv) = 2-Fl^'(-V2, V2;^;il>--->*Jv)|ti=:...=t„=l/fc2, 

whilst the reflection coefficients are given by 
(6.32) 



( ^^N (-V2)iV 2-Fl"'(-V2. iV + 1; ^1, ■ ■ ■ , tiv) 

^! ^F(i)(-i/2,V2;iV;ti,...,t^) 



ti = ...=tjv = l/fe2 



(6.33) 



{N - 1)! lim.^o e2^i''(- V2, -V2; N -l + e;t^,...,tN) 



mi 



F[^\-y2,y2-N-M,---,tN) 



ti = ...=*jv = l/fe2 



In the high temperature phase the diagonal correlation function is 



(6.34) (ao,oaiv,iv) ^ ^^'^A.^f ' k'' ^F'^\y2, 1/2; TV + f ; fi, . . . , t^)!,^^...^,^^^^, 
anci i/ie reflection coefficients are given by 

(-V2)iv 2^'^) (1/2, 3/2; iV + 2; fi, . . . , t^) 



(6.35) 
(6.36) 



(iV + f)!^i.W(i/2,V2;iV + f;ti,...,t^) 



t-\ —... — £ j^ — k^ 



( ,^iv A^! 2J'i''(V2,-V2;jV;ti,...,tiv) 
(V2)iV 27^(1) (1/2, i/2;iV + f;ti,...,tiv) 



Proof. The evaluations in the low temperature phase follow from 14. 51)14. 5111 . although 
some care needs to be taken with fjv because — /i + lj = 0. The limit that arises has a 
series development 



(6.37) lim e^Fl'' (-1/2, - 1/2; iV - 1 + e; ti 

00 

= E 

k:I(,^) = JV 



([-V2ii'')' nf.i(A^ - i + «.) . . . ,fiv) 



[iV] 



(1) 



(iV-1)! 



so that only those terms with lengths l(Hi) = A'^ contribute to the sum. The high tem- 
perature expressions follow from the low temperature ones through the transformation 
/i wi. □ 

It is of interest to note that as A'^ grows more of the leading order terms in the expansion 
of 16.3111 become independent of A'', and the following limit becomes explicit 



(6.38) 



lim (cro.ocrjv.iv) = (1 - A; ^)^/''. 



At zero temperature, k = 00, the solutions simplify to 



(6.39) 



mil 



rN = {N > 1), {cro.ocrjv,]v) = 1, 
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whilst at the critical point, fc = 1, we have the simple solutions 
(6-40) riv= ,7,,,, fN = 



(2iV+ l)(2iV - 1)' ' ' ' 2iV + l' 



AT 



(6.41) {ao,oC7N,N) = Y\ 



r^(i) 



|ir(j + V2)r(j- 1/2)' 

and at infinite temperature they become 

(6.42) = fN = {-!)'' ( Wiv.iv) = (iV > 1), 

in agreement with the known results |42) . 
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